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This paper examines the properties of ‘lattice universes’ wherein point masses are arranged in a 
regular lattice on space-like hypersurfaces; open, flat, and closed universes are considered. The uni¬ 
verses are modelled using the Lindquist-Wheeler (LW) approximation scheme, which approximates 
the space-time in each lattice cell by Schwarzschild geometry. Extending Lindquist and Wheeler’s 
work, we derive cosmological scale factors describing the evolution of all three types of universes, 
and we use these scale factors to show that the universes’ dynamics strongly resemble those of 
Friedmann-Lemaitre-Robertson-Walker (FLRW) universes. In particular, we use the scale factors 
to make more salient the resemblance between Clifton and Ferreira’s Friedmann-like equations for 
the LW models and the actual Friedmann equations of FLRW space-times. Cosmological redshifts 
for such universes are then determined numerically, using a modification of Clifton and Ferreira’s 
approach; the redshifts are found to closely resemble their FLRW counterparts, though with certain 
differences attributable to the ‘lumpiness’ in the underlying matter content. Most notably, the LW 
redshifts can differ from their FLRW counterparts by as much as 30%, even though they increase 
linearly with FLRW redshifts, and they exhibit a non-zero integrated Sachs-Wolfe effect, something 
which would not be possible in matter-dominated FLRW universes without a cosmological constant. 
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I. INTRODUCTION 


Modern cosmology is founded upon the so-called 
Copernican principle , which posits that the universe 
‘looks’ on average to be the same regardless of where one 
is in the universe or in which direction one looks. More 
formally, the Copernican principle states that Cauchy 
surfaces of the universe can be admitted that are homoge¬ 
neous and isotropic, and this symmetry can be expressed 
mathematically by writing the universe’s metric in the 
form 
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where a(t) is the time-dependent conformal scale factor, 
and k is the curvature constant. The sign of k deter¬ 
mines whether 3-spaces of constant t will be open, flat, 
or closed, with k < 0 corresponding to open universes, 
k = 0 to flat universes, and fc > 0 to closed universes; 
one can always re-scale a[t ) and k such that k = +1, 0, 
or —1 as appropriate, in which case a{t) becomes the ra¬ 
dius of curvature for the open and closed universes. This 
family of metrics is known as the Friedmann-Lemaitre- 
Robertson-Walker (FLRW) metric. For the metric to 
satisfy the Einstein field equations of general relativity, 
a(t) must obey the Friedmann equations 
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where A is the cosmological constant, and p and p are the 
energy density and pressure of the fluid filling the space. 
For dust-filled A = 0 universes, the relationship between 
a and t is given by 


a= yp(l — cos ?y) 
t = y(??- sinr/) 

a = (cosh r] — 1) 

t = y(sinh ?7 -r/) 
as depicted in Fig. [Q where 


for k > 0, 


for k = 0, 
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and po is the energy density when a = 1. In the case of 
k > 0, the factor a 0 also corresponds to the maximum of 
a. 

These FLRW models have had great success in explain¬ 
ing much of the universe’s behaviour, including most 
notably the Hubble expansion of the universe, the cos¬ 
mic microwave background (CMB), and baryon acoustic 
oscillations. Indeed, the underlying assumption of ho¬ 
mogeneity and isotropy appears well-supported by pre¬ 
cision measurements showing the CMB to be isotropic 
to within one part in 100,000. Yet recent measurements 
of redshifts from Type la supernovae (SNla), when fit- 
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FIG. 1. Plots of a versus t for open, flat, and closed dust-filled 
universes. 

ted to FLRW models, have led to the conclusion that 
the universe’s expansion is actually accelerating, and to 
account for this acceleration, cosmologists have had to 
modify the standard models by introducing some exotic 
matter, known generally as dark energy. Two of the 
most-favoured dark energy models are an FLRW model 
with non-zero cosmological constant A, which describes 
a negative-pressure fluid acting repulsively under grav¬ 
ity and permeating the entire universe, and a scalar field 
theory called quintessence. However attempts to directly 
observe any exotic matter and determine its nature have 
so far proven unsuccessful, and some cosmologists have 
begun to seek other explanations for the SNla observa¬ 
tions, questioning whether there is any need for exotic 
matter after all. 

There has been an intense debate recently over the rel¬ 
ative importance of inhomogeneities in influencing cos¬ 
mological observations (3-01. In spite of the wide suc¬ 
cess of FLRW models, observations clearly show that the 
late, matter-dominated universe is not homogeneous and 
isotropic except at the coarsest of scales. Instead, matter 
is distributed predominantly in clusters and superclusters 
of galaxies with large voids in between. Standard FLRW 
models ignore such inhomogeneities completely, with pro¬ 
ponents arguing that the inhomogeneities’ effects should 
be small and hence unimportant. Yet others have ar¬ 
gued the converse, proposing that inhomogeneities may 
have non-trivial implications on both the large-scale dy¬ 
namics of the universe and on cosmological observations 
themselves; indeed, inhomogeneities may even explain 
the SNla data as an apparent acceleration without any 
need for exotic matter; this apparent acceleration and 
analogous effects arise as a result of fitting data from an 
inhomogeneous universe onto a homogeneous model. 

Therefore because of the possibly important implica¬ 
tions of inhomogeneities to cosmology, there has been 
much recent effort in trying to quantify their effects. As 
the actual universe’s structure is far too complicated to 
model directly, much effort has been focused on study¬ 
ing the effects of inhomogeneities in toy models. Addi¬ 
tionally, there has been much debate on whether per¬ 
turbed FLRW models can adequately capture the in¬ 
homogeneities’ non-linear structure [J-| 6 (. Thus a non- 
perturbative approach is necessary to unambiguously un¬ 


derstand inhomogeneities. 

Perhaps the simplest inhomogeneous models that can 
account for all known cosmological observations are the 
Lemaitre-Tolman-Bondi (LTB) models [Ml- In such 
models, the observer is located at the centre of a Hubble- 
scale under-dense region in a sphcrically-symmetric, 
dust-dominated universe. These models are therefore 
isotropic but not homogeneous. What is most notewor¬ 
thy is that these models can successfully generate the 
SNla redshift data without requiring any exotic matter. 

FLRW-based inhomogeneous models have also been 
constructed using the ‘Swiss Cheese’ methodology de¬ 
veloped by Einstein and Straus 0]. In these models, 
co-moving spheres are excised from FLRW space-times 
and replaced by either Schwarzschild space-time or LTB 
space-time (21H26|. The resulting models are still exact 
solutions to the Einstein field equations. By construc¬ 
tion, the fluid region surrounding the excised spheres is 
still described by the FLRW metric; therefore these mod¬ 
els would still be dynamically identical to FLRW. How¬ 
ever, there is no a priori reason to believe their optical 
properties should be the same; nevertheless, studies have 
shown that they are generally similar though not identi¬ 
cal to FLRW jaMU. 

In contrast to the above approaches, we shall consider 
a completely different space-time that is not in any way 
based on FLRW nor includes any continuous fluid in its 
matter content. Rather, we shall consider a space-time 
where the matter content consists solely of discrete point 
masses with vacuum in between, as this would provide 
a more realistic description of the late universe’s mat¬ 
ter content. However, for simplicity, our space-time will 
still possess a high degree of symmetry, though not to 
the extent of the FLRW universes. Specifically, we shall 
consider the so-called lattice universes where the mat¬ 
ter content on each Cauchy surface consists of identical 
point masses arranged into a regular lattice. We shall 
focus on lattices constructed by tesscllating a 3-space 
of constant curvature with identical regular polyhedral 
cells; the possible lattices obtainable from such a tessel¬ 
lation are summarised in Appendix [A] To ‘construct’ a 
lattice universe, one of the lattices in Appendix [A] is se¬ 
lected, and all mass in each cell is ‘concentrated’ into 
a point at its centre. The metric consistent with such 
a matter distribution is expected to be invariant under 
the same symmetry transformations that leave the lattice 
invariant, symmetries which include discrete translation 
symmetries, discrete rotational symmetries, and reflec¬ 
tion symmetries at the cell boundaries. In other words, 
the lattice universe should have a metric of the form 

ds 2 = —dt 2 + 7 ^ (t, x) dx a dx b , ( 8 ) 

(o\ 

where 7 ^ b (f, x) is the 3-dimensional metric for constant 
t hypersurfaces, and Latin indices a, b = 1,2,3 denote 
spatial co-ordinates only; the spatial metric 7 ^ (i, x) at 
constant t would possess the lattice symmetries. Effec¬ 
tively, the Copernican symmetries of FLRW universes 
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have been reduced to just these symmetries. Yet these 
lattice universes still satisfy the Copernican principle in a 
coarse-grained manner, in contrast to LTB models where 
homogeneity of the universe is completely surrendered. 

There has been some progress towards determining the 
3-metric 7 ^ (t,x). For the closed universe, the 3-metric 
y( 3 ) (t = 0 , x) on the time-symmetric hypersurface was 
first determined by Wheeler for the 5-cell universe {27] 
and then more recently by Clifton et al. for all closed lat¬ 
tices [28|; this work has been further generalised by Ko- 
rzyriski who examined closed universes with an arbitrary 
number of identical masses, not necessarily arranged in 
a lattice [29|. The time evolution of this initial data has 
been investigated numerically by Bentivegna and Ko- 
rzyriski [30 and analytically by Clifton et al. [31, TO : 
however Clifton et al. only considered the evolution of 
certain highly-symmetric curves on the initial surface, so 
the complete, analytic 4-metric of the closed universe is 
still lacking. The flat lattice universe has been modelled 
perturbatively by Bruneton and Larena [33j and numer¬ 
ically by various authors (Q - l36| : using their perturba¬ 
tive model, Bruneton and Larena have also investigated 
the flat universe’s optical properties [37} • Yet in spite of 
the progress in determining (t, x), a complete, non- 
perturbative solution for any lattice universe remains elu¬ 
sive. Moreover, similar analyses of the open universe, an¬ 
alytical or numerical, have yet to be performed; and the 
optical properties of the closed universe have yet to be 
examined as well. 

Since the complete 4-metric for any lattice universe is 
still unknown, constructing a non-perturbative approxi¬ 
mation of it instead may help provide further insights into 
the universe’s properties. Such an approximation may, 
for instance, offer qualitative insights into the universe’s 
dynamics as well as into the behaviour of photons as they 
propagate along any arbitrary direction; indeed, it would 
be difficult to propagate photons through any universe 
without having a full 4-metric, so for this reason, an an¬ 
alytic approximation of the metric would be especially 
invaluable. In fact, the lattice universe was first studied 
using one such approximation. In 1957, Lindquist and 
Wheeler (LW) [3a [3f| devised a construction wherein 
each polyhedral lattice cell was approximated by a spher¬ 
ical cell with Schwarzschild geometry inside. For this ap¬ 
proximation to work though, the masses were required 
to be of such magnitude and separation that the geom¬ 
etry around each could be reasonably approximated by 
Schwarzschild geometry; for example, should the masses 
get too close together, the deviation around each mass 
from Schwarzschild geometry would become too great 
such that the LW approximation would break down. The 
original construction, however, was restricted to approx¬ 
imating closed universes. Yet even in this limited case, 
Lindquist and Wheeler could obtain a scale factor that 
very closely resembled the scale factor of closed dust-filled 
FLRW universes. Indeed, the LW scale factor would 
asymptotically approach the FLRW scale factor as the 
total number of masses in the universe increased while 


the universe’s total mass was held constant. 

More recently, Clifton and Ferreira (CF) [4(| have re¬ 
visited the LW construction and extended it in several 
notable ways. First, they were able to generalise the con¬ 
struction so that flat and open universes could be mod¬ 
elled as well. Secondly, while Lindquist and Wheeler had 
defined a ‘cosmological time’ co-ordinate in the neigh¬ 
bourhood of the cell boundaries, Clifton and Ferreira 
were able to extend the co-ordinate to cover the entire 
space-time globally. This opened the way to a defini¬ 
tion of co-moving observers for lattice universes, thereby 
allowing quantities such as cosmological redshifts to be 
defined. Finally, Clifton and Ferreira were able to define 
a set of boundary conditions that photons must obey 
when crossing from one cell into the next, thereby al¬ 
lowing photons to be propagated across the entire LW 
model. With their extension, Clifton and Ferreira then 
modelled the flat lattice universe and found that its evo¬ 
lution and redshifts generally agreed well with those of 
the flat dust-filled FLRW universe but that its angular 
diameters and luminosity distances differed [4lj . 

In this paper, we shall further explore the properties 
of the lattice universe using the LW approximation with 
the CF extensions. We shall consider all three types of 
universes, closed, flat, as well as open. First, we shall ex¬ 
tend the LW evolution equation for the closed universe’s 
scale factor and derive analogous expressions for the flat 
and open universes. These scale factors are related to 
the radius of the cell boundary by a constant scaling, 
and in all three cases, the scale factors have the same 
functional form as their FLRW counterparts in Q ([6]). 
Secondly, as Clifton and Ferreira have shown, the LW cell 
boundary radius satisfies a Friedmann-like equation for 
all three types of lattice universes; they have also shown 
that a similar equation holds for the closed universe’s 
scale factor. By using the LW scale factors we derived 
rather than cell radii, we can make the analogy between 
the CF Friedmann equation and the FLRW Friedmann 
equation more salient; specifically, we shall demonstrate 
that the CF Friedmann equation takes a form essen¬ 
tially identical to © for A = 0 and with k in the LW 
case also equalling +1, 0, or —1 according to whether 
the universe is open, flat, or closed respectively. More¬ 
over, we find that the density p in the LW case equals a 
cell’s Schwarzschild mass divided by the spherical volume 
47 rr ^/3 where 77 is the radius of a cell; thus, the density 
behaves as if the cell were a Euclidean sphere. Thirdly, 
the CF global co-ordinate system was found unsuitable 
for studying closed lattice universes. We shall show that 
for closed universes, the CF co-ordinates do not actu¬ 
ally cover the interior of the cell completely but leave a 
gap region instead. For this reason, we must adopt an 
alternative co-ordinate system if we wish to propagate 
photons through closed universes, and in fact, a suitable 
system was proposed by Lindquist and Wheeler them¬ 
selves at the end of their original paper ,38|. Fourthly, 
though our method for propagating photons across the 
universe is strongly influenced by Clifton and Ferreira’s 
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approach, we have adopted several modifications to the 
set of conditions used to propagate photons across cell 
boundaries; these modifications will be justified in the 
paper. Finally, we shall examine the redshifts of pho¬ 
tons travelling along a range of trajectories in closed, 
flat, and open universes. So far, redshifts have only been 
studied in the flat universe using both a perturbative 
approach [37] and the LW formalism [ioi, |4l| ; redshifts 
for the closed and open universes currently do not exist 
anywhere in the literature, so this part of our results is 
entirely new. We shall show that for all three types of 
universes, the redshifts behave broadly in the same way 
as their FLRW counterparts though with certain differ¬ 
ences that can be attributed to the inhomogeneities of 
the lattice. One of our most striking results is that LW 
redshifts can differ from their FLRW counterparts by as 
much as 30%, which has implications on estimating the 
age of the photon’s source, and that an LW universe can 
give rise to a non-zero integrated Sachs-Wolfe effect with¬ 
out needing any cosmological constant. 

This paper is organised as follows. The second sec¬ 
tion will explain the origins of the LW approximation 
and detail its construction; the CF generalisation to flat 
and open universes will also be reviewed. This section 
will also compare the CF and LW global co-ordinate sys¬ 
tem and demonstrate the existence of a gap region if the 
CF system is applied to closed universes. The third sec¬ 
tion will derive the LW analogues of the FLRW scale 
factor and the modified CF Friedmann equation for all 
three types of lattice universes. The latter half of this 
paper will examine the behaviour of redshifts in the lat¬ 
tice universe. The fourth section will briefly explain the 
manner by which redshifts are computed, which follows 
the method of Clifton and Ferreira. The fifth section 
will explain the boundary conditions that are used to 
propagate photons across the boundary between two con¬ 
tiguous cells. The sixth section will explain the numeri¬ 
cal method used to simulate the propagation of photons 
through the lattice. We shall apply Williams and Ellis’ 
Regge calculus formalism for Schwarzschild space-time 
to each cell [42j, [43|; however we have derived an im¬ 
provement to their rules for tracing geodesics through 
the Williams-Ellis Schwarzschild space-time which sig¬ 
nificantly enhances numerical accuracy. Our improve¬ 
ment will be derived in this section with supporting nu¬ 
merical evidence presented in Appendix [B] the improve¬ 
ment is not specific to the models being considered in 
this paper but is applicable to any situation where the 
Williams-Ellis scheme is being used to numerically sim¬ 
ulate geodesics through Schwarzschild space-time. The 
penultimate section presents and analyses the redshift 
results of our simulations. The final section concludes 
with a discussion of possible directions in which this work 
might be extended. 

In this paper, we shall use geometric units where G = 
c = 1. 


II. CONSTRUCTING THE LW 
APPROXIMATION 

The LW approximation of lattice universes was in¬ 
spired by Wigner and Seitz’s [13] method for approxi¬ 
mating electronic wavefunctions in crystal lattices. This 
method approximates the polyhedral cell of a crystal lat¬ 
tice by a sphere of the same volume; any conditions that 
the wavefunction must satisfy on the original cell bound¬ 
ary get imposed on the spherical boundary instead. For 
instance in the original lattice, reflection symmetry at 
the cell boundary means that the wavefunction T of 
the lowest energy free electron must satisfy n ■ V$ = 0 
at the boundary, where the vector n is orthogonal to 
the boundary; in the Wigner-Seitz approximation, this 
same vanishing-derivative condition gets imposed instead 
on the spherical boundary; that is, must now satisfy 
3’F /dr = 0 at the boundary; and this effectively assumes 
the electron potential within a cell to be spherically sym¬ 
metric. The higher the symmetries of the original poly¬ 
hedral cell, the closer the cell resembles a sphere, and the 
more accurate the results obtained from the Wigner-Seitz 
models. Indeed when applied to crystals where exact so¬ 
lutions are known, the Wigner-Seitz construction yields 
very accurate results [H, [if|. 

In analogy with the Wigner-Seitz construction, 
Lindquist and Wheeler approximate each elementary cell 
of the lattice universe by a spherical cell and the metric 
inside each cell by the Schwarzschild metric 

ds 2 = — dt 2 + j^ -^y+r 2 (dQ 2 + sin 2 ddq i 2 ) . 

(9) 

They have called this spherical cell the Schwarzschild- 
cell. Without knowing the true metric, one cannot di¬ 
rectly assess how well is a polyhedral cell in the lat¬ 
tice universe approximated by a spherical one. However 
Lindquist and Wheeler have shown that polyhedra in a 
closed or flat 3-space of constant curvature are reason¬ 
ably approximated by spheres of equal volume [38, 39j ; 
hence when such polyhedra are combined into lattices 
in such 3-spaces, the polyhedral lattice cells would be 
reasonably well-approximated by spherical cells. They 
therefore provide this as partial evidence that the ap¬ 
proximation would probably be reasonable as well for the 
polyhedral cells in the lattice universe. 

In contrast to the Wigner-Seitz lattice, the LW lattice 
is itself dynamical. A test particle sitting at the bound¬ 
ary between two Schwarzschild-cells will by symmetry 
always remain at the boundary. Yet like any other test 
particle in a Schwarzschild geometry, this test particle 
must also be radially falling towards the centre of one of 
the Schwarzschild-cells. And by the same reasoning, this 
particle is also radially falling towards the centre of the 
other Schwarzschild-cell. Therefore this test particle, and 
hence the cell boundary itself, is radially falling towards 
both cell centres simultaneously, as depicted in Fig. [2] 
The boundary’s motion is hence given by the equation of 
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FIG. 2. A test particle co-moving with the boundary between 
two cells must be simultaneously falling towards both central 
masses. Thus the boundary itself must expand and contract 
accordingly. 


FIG. 3. Hypersurfaces of constant Schwarzschild time t inter¬ 
cept each other rather than mesh together at the cell bound¬ 
aries. The normals from different cells point in different di¬ 
rections at the boundary. Therefore, Schwarzschild t cannot 
serve as a global cosmological time co-ordinate for our model. 


motion for a radial time-like geodesic 
2 


2m 


s) =E ~ 1+ — 


( 10 ) 


where r is the proper time of a test particle following 
the geodesic and if is a positive constant of motion; it 
can be shown that \[~E is the particle’s energy per unit 
mass at radial infinity. However this simultaneous free- 
fall of the boundary towards the two masses is actually 
the result of the two masses themselves falling towards 
each other under their mutual gravitational attraction. 
This mutual attraction of all the point-masses thereby 
gives rise to the expansion and contraction of the lattice 
itself, which manifests as the expansion and contraction 
of the Schwarzschild-cell boundary. 

We can see from m that the value of E = Eb at 
the boundary determines whether the cell boundaries will 
expand indefinitely or eventually re-collapse, and hence 
whether the underlying lattice universe is open, flat, or 
closed. For Eb < 1, the boundaries will expand until 
reaching a maximum radius of 


2 TO 



( 11 ) 


before re-collapsing; this corresponds to a closed uni¬ 
verse. For Eb > 1, the boundaries will expand indef¬ 
initely; this corresponds to an open universe. And for 
Eb = 1, the boundaries travel at the escape velocity and 
just reach radial infinity; this corresponds to a flat uni¬ 
verse. Thus it is through this constant Eb that Clifton 
and Ferreira generalise Lindquist and Wheeler’s work to 
flat and open universes in a natural manner. 

In order to ‘glue’ the individual cells together into a lat¬ 
tice, we require that the 3-space of constant time in one 
Schwarzschild-cell mesh at the boundary with the corre¬ 
sponding 3-space of the neighbouring cell. As Lindquist 
and Wheeler have pointed out, two 3-spaces will mesh to¬ 
gether if and only if they intersect their common bound¬ 
ary orthogonally. Surfaces of constant Schwarzschild 
time t do not satisfy this meshing condition; instead, they 



FIG. 4. Hypersurfaces of constant tcf and tlw mesh at 
the boundaries, and the normals from the two cells coincide. 
Therefore either can serve as cosmological time. 


intersect each other when they meet, as illustrated in Fig. 
[3) and some other time co-ordinate must be found. 

Both Lindquist and Wheeler as well as Clifton and 
Ferreira devised new time co-ordinates that satisfied the 
meshing condition. The LW time co-ordinate is defined 
for closed universes only while the CF time co-ordinate 
can only be applied to flat and open universes. Both 
time co-ordinates are defined to equal the proper times 
of a congruence of radial time-like geodesics, which must 
include the geodesics of test particles co-moving with 
the boundary. Test particles following these geodesics 
will travel at the same velocity if at the same radius, 
thus forming freely falling shells. All geodesics’ clocks 
are calibrated to read identical proper time on some ini¬ 
tial space-like hypersurface that intersects the geodesics 
orthogonally. Then if the congruence was well-chosen, 
all other constant proper time hypersurfaces will inter¬ 
sect the congruence orthogonally as well. In particular, 
such hypersurfaces would always intersect the boundary 
orthogonally, thus satisfying the meshing condition, as 
illustrated in Fig. 0] Such a definition of time is well- 
defined globally and can therefore be used as a ‘cosmo¬ 
logical time’ parametrising the lattice universe. Such 
a construction of cosmological time is natural because 
we expect co-moving cosmological observers to be sta¬ 
tionary with respect to constant time hypersurfaces and 
also to follow freely falling geodesics. However there is 
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no unique choice of congruence; any choice that satis¬ 
fies our constraints above is acceptable. Thus herein lies 
the difference between LW time and CF time. Note that 
(fTUl) actually describes any radial time-like geodesic in 
Schwarzschild space-time, including boundary geodesics, 
for which r = r b . Clifton and Ferreira’s congruence con¬ 
sists of geodesics all having the same constant E, namely 
that of the boundary E b . Lindquist and Wheeler’s con¬ 
gruence uses different E < 1 for different geodesics. In¬ 
deed, it may be possible to find an alternative choice of 
congruence based on some stronger physical motivation, 
but we shall not consider this problem here. Instead, we 
shall use LW time for the closed universe and CF time 
for the flat and open universes. 

For boundaries following out-going trajectories, the CF 
time co-ordinate t cf is given by 

jE b -1 + 221 

<1t cf = ^E b dt - ^ _ 2m^ r dr > ( 12 ) 

where E b is the same positive constant as E in (fTUl) for 
the boundary geodesic. The Schwarzschild metric m 
now becomes 


r 



FIG. 5. A qualitative plot of various co-moving test-particle 
trajectories (dashed line) inside a closed Schwarzschild-cell 
and of various constant tcf surfaces (dotted line); the plot 
is of Schwarzschild co-ordinates r versus t. The cell bound¬ 
ary’s trajectory has been drawn with a solid line. The r axis 
passes through the boundary’s moment of maximum expan¬ 
sion, which happens when r CF ~ Tmax■ Since the radial ve¬ 
locity must change sign abruptly when the boundary starts 
closing, there is a noticeable discontinuity in the gradients 
of the test particle graphs at the r axis. The constant tcf 
surfaces for t C f just greater than and less than r maa , are not 
identical. Between the two surfaces is a ‘gap’ region, shaded in 
the figure, where space-time points do not have a well-defined 
cosmological time tcf- 


, 2 1 /, 2m \ 2 2 2 m 

ds = - — 1- dr 2 CF - E b - 1 -|- dr CF dr 

E b \ r ) E b \ r 

dr 2 

+ ^L + r ^ ( m 2 . 

E b 

(13) 

The boundary’s equation of motion is still given by m 
but with the proper time t replaced by t cf , as it can 
be shown that t cf is identical to the boundary’s proper 
time. In CF co-ordinates, the 4-vector tangent to any 
trajectory satisfying (fTUl) for E = E b is 


u a = ^1,^-1 +^,0,0^ . (14) 

For an arbitrary vector n tangent to a constant t cf sur¬ 
face, 


n a = (0, n r , n s , n?) , (15) 

it can be shown that u • n = 0 and hence that surfaces 
of constant t cf are orthogonal to all geodesics satisfying 
m for E = E b , including particularly the boundary 
geodesics. Thus the meshing condition is satisfied, and 
thus t cf can serve as a cosmological time. 

If we attempt to apply the CF co-ordinate system to 
closed universes though, we encounter problems. When 
the boundary is contracting, its radial velocity is given 
by the negative root of (fTUl) instead, and as a result, the 
square-roots in (fl2l) (fT4l) must change sign. However the 
combined co-ordinate patches do not correctly cover the 
interior of the Schwarzschild-cell. Rather, they leave an 
uncovered region centred on the boundary’s moment of 
maximum expansion, as illustrated in Fig. [5] 


We shall denote by T max the boundary’s proper time 
at maximum expansion. For t cf < T max , surfaces of 
constant t cf are generated by the integral curves of m 
for dr CF = 0, that is, by 


dr _ VW( l-^f) 
dt ^E b - 1 + 2s 


\/r b > r > 2m. (16) 


This has a solution of the form t = f(r) + to where to is 
the Schwarzschild time co-ordinate of the surface when 
it passes through r = ro and ro is the Schwarzschild ra¬ 
dial co-ordinate satisfying /(ro) = 0. Thus on the r-t 
plane, like the one shown in Fig. (5j all of these surfaces 
are identical apart from a horizontal shift corresponding 
to different to constants. Increasing to shifts the surface 
to the right. This implies the surface would intercept the 
boundary later in the boundary’s trajectory, and there¬ 
fore t cf would increase as well. As we keep increasing to, 
we shall eventually reach a limiting surface t cf —> T max 
from the right, and this surface we denote by r“ aa; . 

For t cf > T max , surfaces of constant r CF are generated 
instead by the negative of m- This has the solution 
t = —f(r)+to, and again, to is the Schwarzschild time co¬ 
ordinate of the surface when it passes through r = ro . On 
the r-t plane, all r CF > T max surfaces are also identical 
to each other apart from a horizontal shift. Decreasing 
to shifts the surface to the left and decreases t cf . As we 
decrease to, we shall eventually reach a limiting surface 
t C f ~ > T max from the left, and this surface we denote by 

T~*" 

' max 

The complete set of constant r CF surfaces will com¬ 
pletely cover the interior of the cell without any overlaps 
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or gaps if and only if surfaces T max and r^ ax are identical. 
This implies that 

f(r)+t o = ~f{r) + t£ 

=> f{r ) = constant , 

where t// and t// are the to constants for r~ ax and r+ aa ., 
respectively. However we know that function f(r) is not 
constant, so we have a contradiction. 

Instead, we actually have a ‘gap’ between r~ ax 
and T+ ax . Both surfaces meet at {T max , r b {T max )) = 
{Tmax, Tmax)- If we move along T~ ax into the cell, then 
decreasing r would decrease t because dr/dt is always 
positive along this surface. If we move along r+ QX into 
the cell, then decreasing r would increase t because dr/dt 
is always negative along this surface. On the r-t plane, 
T max moves to the left and r+ Q2 , to the right as r decreases 
from r max - Since r~ ax is the rightmost of the t cf < r ma x 
surfaces and r+ aa . the leftmost of the t cf > T max sur¬ 


faces, the region between r“ Qa , and r// ax is not covered 
by any t cf surface and is hence a gap. 

We shall therefore use the LW co-ordinate system, 
which covers closed cells correctly. This system is con¬ 
structed from a congruence of closed geodesics that attain 
maximum radii at the same Schwarzschild time. Being a 
congruence, the geodesics have maximum radii spanning 
the entire range of 2m < r < r b , and from a generalised 
form of CD, the geodesics must therefore have differ¬ 
ent constants E. The geodesics’ clocks are calibrated 
to have identical proper time at maximum radii. The 
LW time co-ordinate t lw at any point is then defined 
to be the proper time of the geodesic passing through 
it. Lindquist and Wheeler also defined a new radial co¬ 
ordinate f LW that is constant along each geodesic and 
equals the geodesic’s maximum radius in Schwarzschild 
co-ordinates. 

The co-ordinate transformation to ( r LW ,r LW ) is then 
given implicitly by 


t = t Q ± 


r LW 1 I / 

2 ^r ) {Tlw 


— r)r 


1 1/2 


2m 


2m In 


■ I r lw 

1 lw — T"0 ^ I TT" 


\ 2m 

r lf2 (f LW /2m - 1) 1/2 + ( f LW - r) 1/2 

r 1/2 (r LW /2m - 1) 1/2 - ( r LW - r) 1/2 
1/2 


( Tlw -A ^ ( Tlw 


Tlw . n \ -l 

- h 2 COS 

2m J \ r L w 


1/2 


2 m 


,,V 2 (~ 

’ \l LW 


\ 1/2 . 
- r) + 


cos 


1 / 2 ' 


(17) 


(18) 


where to and tq are the time co-ordinates at which the 
boundary attains its maximum expansion; the positive 


signs are taken when the boundary is expanding and the 
negative signs when contracting. The Schwarzschild met¬ 
ric now becomes 


ds 2 = —dr: 


Tlw - T 


4 t f lw {f lw / 2,m 1) 


(3 f LW t) 


1/2 


3 Tlw COS 


1 / 2 ' 


-i 2 


df 2 LW +r 2 dfl 2 , (19) 


where r is now a function of t lw and f LW - From this 
metric, it is clear that surfaces of constant t lw are always 
orthogonal to Lindquist and Wheeler’s congruence. Thus 
the meshing condition is also satisfied, and t lw is also 
suitable for use as a cosmological time. 

We shall henceforth drop any subscript labels from the 
cosmological time r. In the case of the closed universe, 
r will denote t lw , while in the flat or open universe, it 
will denote t cf . 


III. THE COSMOLOGICAL SCALE FACTOR 
AND THE FRIEDMANN EQUATIONS 

There is a clear analogy between the size of the 
Schwarzschild-cell 7>,(r) and the FLRW scale factor a(t), 
as both provide a scale of their respective universe’s size. 
Indeed we expect the lattice universe’s scale factor a(r) 
should be a function of rt, alone, that is a(r) = a(rb(r))o 


1 We shall use a(t) when we refer to the FLRW scale factor and 
a(r) when referring to the lattice universe scale factor. The two 
functions are completely independent. 
















If we expand rb, then each cell of the lattice will expand 
by some scale £, and therefore the lattice as a whole will 
expand by £. We shall take a(r) to be related linearly 
to r’i(r), that is, a(rb{r)j = ar&fY) for some constant 
a > 0. As we shall see, a(r) then depends on r in a 
manner analogous to how a(t) depends on t for FLRW 
universes. 
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FIG. 6. Rescaling the cell size by £ rescales the entire lattice 
universe by £ as well. 


Lindquist and Wheeler were able to derive a paramet¬ 
ric relationship between a(r) and r closely resembling 
that of GD, thus showing that a(r) behaves essentially 
identically to its FLRW counterpart a(f); the only dif¬ 
ference was the factor of ao- However, their results were 
limited to the closed lattice universe. We shall generalise 
Lindquist and Wheeler’s results, showing that this iden¬ 
tical behaviour between a(r) and a(t) holds for lattices of 
all curvatures. Equation m can be integrated to obtain 

t = (V r (r E - r) +r E cos -1 ^J ~-- j ( 20 ) 

for 0 < E < 1, 



for E> 1, 


where 



for E ^ 1 and E > 0, 
for E= 1. 


(23) 


In the case of E < 1, r E corresponds to the maximum 
radius attained by the geodesic, and for the boundary 
geodesic, where E = E b , this maximum radius is identical 
to m • For boundary geodesics, where r = r b , equations 


(l20l) to (l22l) can be re-cast into the form 


r b = —~ (1 - cost?) 


r b = 



for 0 < E b < 1, (24) 


for Eb = 1 , 


(25) 


r b = (cosh 77 - 1 ) 


r Eb r Eb 

T = ^“V^ (sinhr?_7?) 


for Eb > 1 , 


(26) 


where r Eb simply denotes r E when E = E b and 77 is simply 
a parametrisation. If we choose a to be 


(27) 


a = (1 — cos 77 ) for 0 < E b < 1, (28) 

2 V 2 m 



2m’ 

then the scale factor a(r) is given by 

7- 

1/3 


a = tV3 f ° r Eb < (29) 


a = —— i —— (cosh 77 — 1) for E b > 1. (30) 

2 V 2m 


Comparing the relations just obtained for a(r) and r with 
their counterparts in (U) ([ 2 ]) , we find that a(r) and r have 
the same functional form as their FLRW counterparts for 
all background curvatures. Moreover, we have deduced 
an expression for the factor a that is the same for all 
universes. The only difference between the lattice and 
FLRW relations is that oq for the lattice universe is 



Equivalently, by equating (l3lT) with J7|), we can say that 
the density po f° r the lattice universe is 


m 



(32) 


The denominator is simply the ‘Euclidean volume’ of a 
Schwarzschild-cell with radius r b corresponding to a = 1. 

As Clifton and Ferreira have noted, the radius of the 
cell boundary r&(r) satisfies an equation strongly resem¬ 
bling the Friedmann equation © for A = 0; they also 
found a similar relation for the closed LW universe’s scale 
factor. By consistently using the LW scale factors just de¬ 
rived rather than t-j,(t), we can re-express the CF Fried¬ 
mann equation in a form that makes its resemblance to 
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([2]) much more salient; the CF Friedmann equation then 
becomes 


a(r )\ 2 87 xp k 

a(r )) 3 a(r) 2 ’ 


(33) 


where N is the total number of cells in the lattice. For 
hyperbolic spaces, we define if; analogously in terms of hy¬ 
perbolic angles. If ro is the radius of one such spherical 
cell, then we can relate tq to the comparison hypersur¬ 
face’s radius of curvature a LW by 


where k = a 2 (l — Eb) plays the role of the curvature 
constant for the lattice universe, and density p is given 

by 


- _ 171 _ 171 /o 4 \ 

P (| 7 ra(r) 3 )/a 3 |tt rf ‘ 

As noted earlier, Eb determines whether the universe will 
be open, flat, or closed, and we see that k will take on 
the correct sign accordingly. In fact, if we make use of 
( 1751 ) and ( 1271 ) . k simplifies to — 1 ,0 , +1 for open, flat, and 
closed lattice universes, respectively, much like its FLRW 
analogue. 

To complete our model, all that remains is to specify 
a or Eb. However, Lindquist and Wheeler have actually 
provided an argument to specify a and an independent 
argument to effectiveljH specify £)What is intriguing 
is that their choices of a and Eb are consistent with the 
relationship between the two quantities required by ( 1771 ) . 
We shall summarise Lindquist and Wheeler’s arguments 
leading to their choices. Although they dealt only with 
closed universes, we shall, in this paper, generalise their 
arguments to flat and open universes as well. 

We begin with Lindquist and Wheeler’s choice of a. 
Note that hypersurfaces of constant t in the FLRW met¬ 
ric © correspond to 3-spaces of constant curvature, and 
through an appropriate choice of scaling, the radius of 
curvature is given by the scale factor a(t). We embed 
the same type of lattice as the lattice universe on such a 
hypersurface with the appropriate curvature @ We shall 
call this hypersurface the comparison hypersurface. We 
then approximate each polyhedral cell by a sphere of the 
same volume 0 For closed hyperspheres, we define 0 to 
be the angle between the centre of the spherical cell and 
its boundary as measured from the centre of the hyper¬ 
sphere. Then 0 is given implicitly by 


1 20 — sin 20 

N 2t r 


2 Recall that the quantity E 5 was introduced to the LW construc¬ 
tion afterwards by Clifton and Ferreira. However, Lindquist and 
Wheeler imposed a tangency condition on the Schwarzschild- 
cells, which we shall soon describe, that effectively specifies E 5 . 

3 That is, for closed lattice universes, we use hyperspheres; for flat 
lattice universes, Euclidean space; and for open lattice universes, 
hyperbolic space. 

4 This is actually Clifton and Ferreira’s generalisation. Lindquist 
and Wheeler’s original condition was that the spheres occupy 
N —1 of the comparison hypersphere’s total volume, where N is 
the total number of cells. Such a condition clearly needs to be 
modified for flat and open universes, as both N and the hyper¬ 
surface volume are infinite in these cases. 


a L w 


^0 

x(V0 ’ 


where the function x(V 0 is given bjfi 


xW 


{ sin'i j) for closed universes, 
1 for flat universes, 

sinhi/> for open universes. 


(36) 


(37) 


Lindquist and Wheeler identified the Schwarzschild-cell 
radius of the lattice universe with ro and used the 
corresponding a LW as given by ( 1351 ) to be the lattice uni¬ 
verse’s scale factor; that is, 


i \ _ r &( T ) 

aLw(T) ~ m- 

Hence they have chosen a = 1/xWO- 
Eb would therefore correspond to 


(38) 

According to ( 1771 ) . 


cos 2 1 /> 

for closed universes. 

1 

for flat universes, 

cosh 2 if) 

for open universes. 


(39) 


However, Lindquist and Wheeler prescribed Eb in¬ 
dependently of m- They instead embedded the 
Schwarzschild-cell in the comparison hypersurface and 
required it to be in some sense tangent to the hyper¬ 
surface. Their prescription for the embedding is as fol¬ 
lows. Suppose we replaced each of the original spherical 
cells in the hypersurface with a Schwarzschild-cell from 
the lattice universe. Then we want to choose an Eb to 
make the Schwarzschild-cell tangent to the hypersurface. 
Lindquist and Wheeler formulate this tangency condi¬ 
tion as follows. Take any great circle on the boundary of 
the original sphere and compare its circumference with 
that of the corresponding great circle on an infinitesi¬ 
mally smaller sphere. Depending on which hypersurface 
the cell is embedded in, the circumferences will obey the 
relation 

1 d(circumference) 

27t d(radial distance) 


5 Note that when k = 0 in the FLRW metric 0 . we have complete 
freedom to make a re-scaling of the form r —> and a —»• a /£. 

Hence we can always choose a comparison hypersurface such that 
cllw = tq- 
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1 d(2ira LW sin^) 
2n a LW di,b 

< 1 

1 d(27ra iw sinh'0) 
„ 27r a LW dtp 

( COS Ip, 


for closed hyperspheres, 
for flat hypersurfaces, 
for open hypersurfaces, 


r 

Lcoshi/), 


as depicted by Fig. [7] to Fig. [9] Then for the 
Schwarzschild-cell to be tangent to the hypersphere, we 
also require that at its boundary, 


cos ip for closed universes, 


1 d(circumference) , . _ 

--- P- = < 1 for fiat universes, 

27r d(radial distance) 

L cosh 7/; for open universes. 

(40) 

From the Schwarzschild metric expressed in CF co¬ 
ordinates m we can see that d(circumference) = 2-Kdr b 
and that d(radial distance) = dr b /y/Ep- Thus, 


1 d(circumference) 
27 t d(radial distance) 


yfW b . 


From the Schwarzschild metric in LW co-ordinates m, 
it can also be shown, by making use of m, that 
(27r) _1 d(circumference)/d(radial distance) is the same. 
Solving for E b , we then obtain the same result as in (l39l) . 

We close this section with a few remarks about the 
large-./V behaviour of the closed LW model. By taking 
N —^ oo, we can deduce what limiting behaviour the 
model would approach as the number of masses increases 
while the universe’s total mass is held constant. It can be 
shown from (l35l) that as N —> oo, the angle xp —> ■ 

If M = Nm is the total mass of the universe, then 


o M 2 M 

lim a m = lim —-=— =-, 

N—>oo N—too N sin ip 37 t 

and therefore the density po for the lattice universe, as 
given by (1321) . becomes 


lim po 

N—t oo 


M 

2n 2 


(41) 


In closed FLRW space-time, a hypersphere of constant 
t has a volume of 2n 2 a(t) 3 . Since the FLRW po is de¬ 
fined to be the density when a(t) = 1, then po also equals 
M/27T 2 . Therefore as N —> oo, the lattice universe den¬ 
sity po approaches its FLRW equivalent p 0 , and hence the 
lattice universe factor a o approaches its FLRW equiva¬ 
lent as well. Consequently, a(r) in (1751) approaches a(t) 
in ([4]), and r in (l24l) becomes identical to t in d4|. We 



FIG. 7. A spherical cell of radius aLwsimj) and an infinites¬ 
imally smaller shell, indicated by the dashed line, have been 
embedded in a three-dimensional hypersphere of radius a LW , 
with one of the angular dimensions suppressed. The radial 
distance between the two shells, as measured along the 3- 
sphere, is a,Lwdil>. Because of the curvature of the underlying 
3-sphere, we have that d(circumference)/d(radial distance) = 
27T cos ip. 



FIG. 8. An analogous embedding in a flat hypersur¬ 
face. In this case, the circumference of the boundary is 
27r a,Lw', the radial distance between the two shells is dr; and 
d(circumference)/d(radial distance) = 277. 



FIG. 9. An analogous embedding in a hyperbolic hypersur¬ 
face. The circumference of the boundary is 2-iraLw sinh ip; 
the radial distance between the two shells is a,Lwdip; and 
d(circumference)/d(radial distance) = 277 cosh ip. 


also note that p in (l34l) becomes 


lim p 
N—*oo 


M 

2i r 2 a 3 ’ 


(42) 


which is identical to the FLRW energy density p. This 
implies that the CF Friedmann equation (l33l) would be 
identical to the FLRW Friedmann equation @ for k = 1 
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and A = 0. Therefore, as N increases and the closed 
lattice universe approaches the continuum limit, it be¬ 
comes increasingly similar to the closed dust-filled FLRW 
universe, which is consistent with the matter content it¬ 
self becoming increasingly similar to homogeneous and 
isotropic dust. In fact, Korzynski has proven that this 
large -N correspondence between the two closed universes 
is indeed true for the exact solution of the Einstein Held 
equations, at least on the hypersurface of time-symmetry 
[29| . He actually considered universes where there could 
be an arbitrary number N of identical Schwarzschild-like 
black holes and demonstrated that as N —> oo, the hyper¬ 
surface approaches its FLRW counterpart exactly, both 
in terms of the scale-factor and in terms of the energy 
density. Therefore in this regard, the LW approximation 
agrees with the exact solution. 

IV. REDSHIFTS IN THE LATTICE 
UNIVERSE 

Cosmological redshifts, 1 + z FLRW , in FLRW universes 
are defined with respect to sources and observers that are 
co-moving with the universe’s expansion. The redshifts 
are given by 

1 + z FLRW = —, (43) 

a e 

where a 0 and a e are the cosmological scale factors at the 
moments of observation and emission, respectively. 

We should like to find the lattice universe analogue to 
Zflrw so that we can compare redshifts in the two types 
of universes. In general, redshifts 1 + z are defined as the 
ratio of the photon frequency measured at the source to 
the frequency measured by the end observer. The clos¬ 
est analogy to co-moving sources and observers in the 
lattice universe would be sources and observers that are 
co-moving with respect to a constant r surface. All such 
observers would be following radial geodesics that obey 
(Hhd, with E fixed to be Ef, in the flat and open universes 
but geodesic-dependent in the closed universe. If u is 
the 4-velocity of an observer and k the 4-momentum of 
a photon as it passes the observer, then the photon fre¬ 
quency measured by the observer would be —u • k. For 
the lattice universe, the cosmological redshift is therefore 
given by 


where the subscripts s and o denote ‘source’ and ‘ob¬ 
server’, respectively. Clifton and Ferreira have con¬ 
strained their consideration to observers that are co¬ 
moving with the photon’s source; that is, at any time r, 
the observer would be at the same radius as the source in 
their respective cells’ Schwarzschild co-ordinates; this has 
been illustrated in Fig. [TQ1 To facilitate comparison with 
Clifton and Ferreira’s results, we shall compute redshifts 



FIG. 10. A photon travels from radius Ti at cosmological time 
n in one cell to radius r/ at time r/ in the next, as indicated 
by the long solid arrow. We assume that the photon always 
passes through the boundary in the manner illustrated here: 
that is, the point of crossing is a point of tangency between the 
boundaries of the two cells. A photon is emitted at (n, ri ) by 
a source travelling along a geodesic given by EB and observed 
by an observer in another cell. Although in a different cell 
from the source, the observer is still ‘co-moving’ with the 
source; that is, for any r, the observer is always at the same 
radius as the source in their respective cells, and this requires 
the observer to travel along a geodesic with the same E as 
the source. 

for the same set of observers. 

Following Clifton and Ferreira’s example, we shall also 
use the lattice universe’s scale factor a(r) at the moments 
of emission and observation to compute z FLRW for com¬ 
parison; we thus re-express (l43l) as 

l + z FLRW = ^\, (45) 

rb{T e ) 

where we have made use of the fact that a(r) = ar^r). 

V. PROPAGATING PHOTONS ACROSS 
CELL BOUNDARIES 

In order to propagate photons through the lattice uni¬ 
verse, we must first specify what boundary conditions 
trajectories must satisfy whenever they pass from one 
cell into the next. Before discussing the boundary con¬ 
ditions though, we first note a difference between Clifton 
and Ferreira’s choice of boundary geometry and ours. 
Clifton and Ferreira converted from spherical cell bound¬ 
aries back to polyhedral ones, deducing the polyhedral 
boundary velocity from the requirement that the spher¬ 
ical cell always have the same ‘Euclidean volume’ as the 
polyhedral cell. We shall instead continue to use spheri¬ 
cal boundaries and propagate photons across boundaries 
in the manner illustrated in Fig. 1101 that is, wherever 
a photon crosses, we always regard the point of cross¬ 
ing as a point of tangency between the two neighbour¬ 
ing cell boundaries. We assume that the boundaries are 
tangent not just in (3+l)-dimensional space-time as a 
whole, but also in each 3-dimensional constant-r hyper- 
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surface, which would be the case had we been using the 
true polyhedral lattice cells. As Lindquist and Wheeler 
have argued, spherical boundaries should be a good ap¬ 
proximation to the shape of polyhedral boundaries, with 
the approximation improving as the number of symme¬ 
tries increases. Therefore any errors due to this approx¬ 
imation would be small to begin with. We further ar¬ 
gue that although spherical cells may tile the lattice uni¬ 
verse with gaps and overlaps, an arbitrary photon would 
on average travel through an equal number of gaps and 
overlaps such that the overall error approximately can¬ 
cels out, as illustrated in Fig. [XT] Additionally, we note 



FIG. 11. A flat lattice universe with cubic cells approximated 
by spheres. The central masses are depicted by dots. A pho¬ 
ton will on average travel through an equal number of gap 
and overlap regions. Because we are using spherical bound¬ 
aries, at boundary crossing A, the photon will jump over a 
gap region, but at crossing B, it will pass through the over¬ 
lap twice. The cubic boundaries are tangent to each other 
everywhere along the boundaries while the spherical bound¬ 
aries are nowhere tangent to each other except where they 
intersect; though even there, the tangency is only in the 
entire (3+l)-dimensional space-time, not in any specific 3- 
dimensional constant-r hypersurface. However in the LW 
approximation, we shall assume that the spherical bound¬ 
aries are everywhere tangent to each other, both in the en¬ 
tire (3+l)-dimensional space-time and in each 3-dimensional 
constant-r hypersurface; this is because we are treating the 
spherical boundaries as if they were approximately identical 
to the cubic ones. 

that in Clifton and Ferreira’s polyhedral cells, the mesh¬ 
ing of constant-r hypersurfaces at cell boundaries is lost, 
with the hypersurfaces now meshing in an average man¬ 
ner only. In contrast, the hypersurfaces in the original 
spherical cells, by construction, mesh exactly. Finally, 
the idea of replacing polyhedral boundaries completely 
with a spherical approximation seems closer in spirit to 


the original idea of Wigner and Seitz. Therefore accord¬ 
ing to our choice of boundary geometry, if a photon exits 
its current cell at a radius ri = r^, we require it to en¬ 
ter the next cell at radius r 2 = rb . Because of cell 2’s 
spherical symmetry, we are free to choose any 9 and <f> 
co-ordinate for the entry point. 

We now present the conditions that photon trajectories 
must satisfy when propagated across boundaries. These 
conditions are applied locally at any pair of exit and en¬ 
try points. Our conditions are founded upon Clifton and 
Ferreira’s principle that any physical quantity should be 
independent of which cell’s co-ordinate system an ob¬ 
server co-moving with the boundary may choose to use. 
In the context of photon trajectories, we require that 

1 ) the photon frequency match across the boundary, 

ui ■ ki = u 2 • k 2 , (46) 

2) and the projection of the photon’s 4-momentum 
onto the vector n r orthogonal to the boundary 
match across the boundaryjj 

■ ki =-n£ • k 2 ; (47) 

it can be shown that ( n r ) a = ( 0 , \/Eb, 0 , 0 ) for 
both the LW and CF co-ordinate systems. There 
is a negative sign in the above equation because 
n r always points radially out of its respective cell, 
whereas the radial direction of k would be out of 
one cell and into the next. 

These conditions along with the normalisation k • k = 0 
are sufficient to deduce the components of k 2 in terms of 
ui, u 2 , and ki. 

The vectors u and n r imply a decomposition of k into 
the form 

k = —(k • u)u + (k • n’’)n r + k 22 , (48) 

where it can be shown that k 22 ■ u = k 22 ■ n r = 0. We 
note that because of the cell’s spherical symmetry, we can 
always choose polar co-ordinates such that k lies in the 
9 = 7 t/2 plane, thereby allowing us to suppress the 9 co¬ 
ordinate. In this case, k 22 would take the form ( k n ) a = 
(0,0, 0,fc n ). The above conditions imply that (fc ) 2 ) 2 = 
(/c^ 2 ) 2 , and by a suitable choice of the <p co-ordinates, we 
can always make fcp = fc^ 2 . We also note that there is no 
physical reason for the photon trajectory to refract when 
passing through a boundary, and therefore the 9 = n/2 
planes of the two cells should be aligned. 

Using the above decomposition and the two bound¬ 
ary conditions, we can therefore express the photon’s 4- 


Because of the normalisation condition, this second requirement 
is equivalent to requiring that the space-like projection of k onto 
the boundary be identical in both cells. 
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momentum k 2 in the new cell as 

k 2 = -(ki • m)u 2 - (ki • n[)nz + k^. (49) 

VI. NUMERICAL IMPLEMENTATION OF 
THE LW MODEL 

To numerically simulate photons propagating through 
the lattice universe, we have chosen to implement 
Williams and Ellis’ Regge calculus scheme j42| for 
Schwarzschild space-times. We shall use their method to 
discretise the space-time of each Schwarzschild-cell into 
Regge calculus blocks, and then employ their geodesic¬ 
tracing method to propagate photons and test particles 
through a cell, though with some modifications. In gen¬ 
eral, the discrete space-times of Regge calculus are ex¬ 
pected to converge at second order in the lattice edge- 
lengths to the continuum space-times of general relativity 

m 

Under Willliams and Ellis’ scheme, Schwarzschild 
space-time is discretised as follows. A grid is constructed 
in the Schwarzschild space-time such that along any 
particular gridline, only one Schwarzschild co-ordinate 
(t, r, 0, <fi) changes while the other three are held con¬ 
stant. The lines intersect at constant At, A r, A8, or A 4> 
intervals, thus forming the edges of curved rectangular 
blocks. Each of these blocks gets mapped to flat rect¬ 
angular Regge blocks such that the straight edges of the 
Regge blocks have the same lengths as the curved edges 
of the original blocks; Figure [T5] illustrates an example 
of a curved Schwarzschild block with its corresponding 
Regge block. The Schwarzschild co-ordinates of the orig¬ 
inal blocks’ vertices (t,, rj, 6k , (f>i ) are now taken over as 
labels for the Regge blocks’ vertices. As the Regge blocks 
are flat, the metric inside is simply the Minkowski metric. 

Suppose a Regge block had vertices (ti, r,j , Ok , 4>i ) 
and (ti+i,rj+i,0fc+i,0j+i); then we shall use the label 
(ti, rj,8k, 4>i ) to refer to this block as a whole. The edge- 
lengths for this block are 


d[(ti,rj,9 k ,<j>i), (t i+ i,rj,8 k , (jn)\ 


d[(ti,r 3 ,8 k ,(t>i), (ti, rj + i,9 k , <f>i)] 



d[(ti,rj,9 k ,<t>i), ( ti,rj,8 k+ i,(t>i )] = rj(8 k +i - 8 k ), 


d [(ti, rj , 8 k , (fn) , (U , rj , 9 k , (j>i +1 )] = rj sin 8 k (<f>i +1 - cjn). 

Since photons always move in planes, we can always 
suppress one of the angular co-ordinates by setting 9 to 
be 7 t/ 2. We then set-up a Minkowski co-ordinate sys¬ 
tem (t,p,i/j) in the block, with block vertices located 
at (±Ti,-pi,±%pi) and (± r i+1 , pi, ± V^+i), where n, p u 


and *l>i are defined b£| 


(50) 

2-1 1/2 


and wher<0 



Williams and Ellis propagate geodesics through Regge 
blocks on the principle that geodesics should follow 
straight lines both within a block and on crossing from 
one block into the next. There is an apparent refraction 
of the geodesic in crossing into a new block because the 
(r, p, %ji) co-ordinate systems of the two blocks are not 
aligned; therefore the same tangent vector of a geodesic 
would be represented differently in different blocks’ co¬ 
ordinate systems. Williams and Ellis have demonstrated 
that their scheme successfully reproduces the orbits and 
redshifts of particles travelling in Schwarzschild space- 
time Big. 

Referring to Fig. 1121 we now summarise the rules for 
propagating geodesics from one block into the next. Sup¬ 
pose the particle is at position (to, po, ipo) and travel¬ 
ling in direction k a = (k\,k 2 ,ko). It exits the block at 
(t',p',i/j') = (To,po,ipo) + Ak, where the value of A de¬ 
pends on the face exited. 

1) If the particle exits the block by the top face, then 

> _ t - t 0 + po tanh & 
k\ — k 2 tanh /3j ’ 

where t = (r^+i + Ti )/2 and tanh/% = (rj+i — 
Ti)/2pi. In the new block, the particle’s new tra¬ 
jectory is given by applying to k the matrix 

/ cosh2/Ij — sinh 2/3,; 0\ 

I — sinh 2Bi cosh 2/3, 0 ] , (53) 

V 0 0 ' V 

and the particle’s starting position is given by 



2*02 — ^ i (0J+1 0^)5 

_’j 2 , (n+i - n) 2 (ipi+i - i/i) 

Pi 1 , , 


7 Note that an arbitrary index j appears in the expressions for Ti 
and However, because both (tj+\ — tj) and (0j+i — <j>j) are 
actually independent of j , it does not matter what value of j is 
used to calculate Ti and 

8 In O , there was a missing factor of 2m in front of the logarithm 
for this equation; this factor has been restored here. 
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FIG. 12. On the left is an example of the original Schwarzschild block with one angular co-ordinate suppressed, and on the 
right is the Regge block to which it is mapped, (r, p, i/j) is a Minkowski co-ordinate system within the Regge block. 




2) If the particle exits by the back face, then 


A = 


ip — V’o + Po tan a, 

k 3 — k 2 tan a* 


(54) 


where ip = (ipi +1 + ipi)/2 and tana, = (t/>i+i — 
tpi)/2pi. In the new block, the particle’s trajectory 
k gets transformed by 

/f ° 0 \ 

0 cos2ai sin 20 ^ I , (55) 

yO — sin 2 on cos 2 on) 

and the particle’s starting position is (t',p', — ip'). 

3) If the particle exits by the front face, therp] 


A = 


-ip - -00 ~ potanaj 
&3 + k 2 tan on 


(56) 


where ip and tanaj are the same as above. In the 
new block, the particle’s trajectory k gets trans¬ 
formed by the same rotation matrix as (l55l) but 
with angle — ct^ instead. The particle’s starting po¬ 
sition in the new block is again (t ', p' ,—ip r ), since 
ip' < 0 at the front face of the original block. 

4) If the particle exits by the right/left face, then 


_ ±Pi ~ Po 

k 2 


(57) 


and p' is simply p' = ±pi . There is no refraction as 
the particle enters the next block. Its new starting 
position is (r', —pi + i, ip') if entering the block to 
the right and (r', Pi-\,ip') if entering the block to 
the left. 


9 In [HI, there were a few sign errors for this equation, which have 
been corrected here. 


However from numerical simulations, we have discov¬ 
ered an empirical relation for the tangent 4-vectors of 
radially out-going time-like geodesics. This has led us 
to introduce a correction to the above propagation rules. 
We have found that as a function of r t . these 4-vectors 
obey 

Ke 99 e = 



(58) 

where f mQa; is a constant of motion. Our numerical re¬ 
sults supporting this have been presented in Appendix 

E 

This relation can be derived analytically by compar¬ 
ing u Kegge with its counterpart 4-vector in continuum 
Schwarzschild space-time. The tangent 4-vector of a radi¬ 
ally out-going time-like geodesic in Schwarzschild space- 
time is 

= (i, r, sV) 



where r max is also a constant of motion and is related to 
E in (flOll by r ma x = 2to/(1 — E). For in-going geodesics, 
the radial components of both u Hegge and u Sch „ arz , would 
have an additional negative sign. If we take the scalar 
product of u flegge and u Schwarz with the unit vector in 
the time direction, that is, with f a — (1,0,0) for Regge 

space-time and t a = ^(1 — 2m/r)~ 1 ^ 2 , 0 , 0 ^ for contin¬ 
uum Schwarzschild space-time, we have that 
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and that 


of motion E , which we define by the relation 


t 


u 


Sch- 





If we identify r mQX with fmax, then these two expres¬ 
sions are identical whenever r = n in the continuum 
space-time. Similarly, if we take the scalar product 
with the unit vector in the radial direction, that is, 
with p a = (0,1,0) for Regge space-time and f a = 

^0, (1 — 2m/r) 1 ^ 2 ,0^ for continuum space-time, we have 
that 


P ' H-Regg< 



and that 


( 2m _ 2 m 

Vi Vmax 

l 2m 

Again if we identify r max with r max , then the two expres¬ 
sions are also identical whenever r = in the continuum 
space-time. Therefore provided f max = rmax-, we see that 
(1551) is indeed the Regge analogue of (1551) . Furthermore, 
the choice of r mQX and r max determines whether the re¬ 
sulting particle orbit will be closed or open in the corre¬ 
sponding space-time. If 2 m/r max > 0, then the orbit in 
the continuum space-time will be closed and r max would 
indeed be the maximum radius of the orbit. Similarly, 
we found that if 2 m/f m ax > 0, then the orbit in Regge 
space-time will also be closed, and the maximum radius 
would be fmax■ If ‘Zm/rmax = 0, then the geodesic will 
just reach spatial infinity in the continuum space-time, 
and (l59l) gives the particle’s escape velocity as a func¬ 
tion of its radial position r > 2m. Similarly, we show 
in Appendix iBl that when 2m/fmax = 0, then (1551) gives 
the escape velocity for a test particle in the Regge space- 
time. Finally if 2 m/r m ax < 0, then the orbit in con¬ 
tinuum space-time will be open. And if 2m/fmax < 0, 
then the orbit in Regge space-time will also be open, as 
we show in Appendix [Bj We shall henceforth make the 
identification of r mQX = r max - 

Inspired by this, we can generalise the expression for 
u Regge to any geodesic in Regge Schwarzschild space- 
time. The Lagrangian for particles moving in continuum 
Schwarzschild space-time can be written as 



L=- 




r 2 Q 2 , 


where the dot denotes differentiation with respect to 
some parameter A. Since 0 = dL/dt , we have a constant 


E := 



(60) 


Similarly, since 0 = dL/dfl, we have another constant of 
motion J defined by 


J := r 2 n. (61) 

These two constants correspond to the square of the par¬ 
ticle’s energy per unit mass at radial infinity and to the 
angular momentum. For radial time-like geodesics, E 
here is the same as E in (fTOl) . In terms of these con¬ 
stants, the tangent 4-vector vschwarz to the particle’s 
geodesic can be expressed as 



which is clearly a function of r alone, since f can be 
deduced from t and ff through the normalisation of 

^Schwarz • 

Let v Regge denote the Regge analogue of v Schwarz . If we 
assume that v Regge and v Schwarz are related in a manner 
analogous to how u Regge and u Schwarz are related, then 
we can deduce the components of 'v Regge from v Sch „ ars ,. 
Specifically, by equating f-v Regge to t-u Schwarz , we deduce 
t to be 


Ve 

T ~ (i -2m) 1/2 ' 

Similarly by equating /p ■ v Regge to fl ■ u Schwarz , where 
ip a = (0, 0, 1) and fi a = (0, 0, r -1 ), we deduce ip to be 



Thus the components of v Re99e are given by 

= 7 ) ’ ^ 

with p deducible from the normalisation of v Re39e . Fi¬ 
nally, it is easy to verify that p ■ v Regge and r • v Schwarz 
are then consistent. 

Our generalised expression for tangent 4-vectors can be 
used for both null and time-like geodesics following both 
radial and non-radial trajectories. For time-like radial 
geodesics, where J = 0, it can be shown that v Hesge re¬ 
duces to n RBgge . As we wish to simulate geodesics in con¬ 
tinuum rather than Regge Schwarzschild-cells, we have 
introduced a correction to the Williams-Ellis scheme: 
whenever a geodesic crosses a right or left face, its tan¬ 
gent 4-vector is changed back to HMD- 

Having established the modified Williams-Ellis scheme, 
we now discuss its application to the Schwarzschild-cell. 
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The boundary is simulated by propagating a test parti¬ 
cle that is co-moving with the boundary. Whenever the 
photon reaches the same radial and time co-ordinates as 
the boundary particle, then at that point the photon ex¬ 
its one cell and enters the next, and we need to apply 
conditions (|46l) and (S3 to determine the photon’s new 
trajectory. However vector components will differ from 
previously as we are now working in Regge space-time. 

Let u a = (u T ,u p , 0) be the boundary particle velocity 
and k a = (fc r , k p , k^) be the photon’s 4-momentum at the 
boundary. We shall use subscripts 1 and 2 to indicate the 
cell being exited and the cell being entered, respectively. 
As in (l48l) . we decompose the Regge vector k into 

k = i'll + n p + n^, (64) 

where v = u • k, ( ri^) a = (0,0, ri^), and n p satisfies 
n p • u = n p • = 0. If n p is the normalised vector of n p , 

that is n p = n n p , then the components of n p are 

(n p ) a = (w p , u T , 0), 

as this satisfies all orthogonality relations required of n p . 
We then deduce n to be 

n = n p • k = ( u T k p — u p k T ). (65) 

Condition (1461) implies that 

Vi = v 2 - ( 66 ) 

Condition (1471) in this context is equivalent to 

k| • n p = —k 2 • n p 2 , 

which implies that 

ni = —ri2- (67) 

As at the end of Section El these conditions imply that 

(nf) 2 = (n 2 ) 2 , and we again have freedom to choose 
polar co-ordinates such that 

nf = n$, (68) 

which we shall henceforth assume to be the case. Using 

relation (l65l) . conditions (l66l) and (l67l) . relation (1681) . and 
decomposition (|64)l . we can express the components of 
k 2 as 

k% = (vu2+(uikl~u\ki)u2, vu^+iuikl—ulki)^, kf). 

(69) 

We note that constant E in (17>3l) will differ between ki 
and k 2 , but J will remain the same since k^ = J/r is 
identical on both sides of the boundary. 


VII. REDSHIFTS FROM LATTICE 
UNIVERSE SIMULATIONS 

We have simulated the propagation of photons through 
multiple Schwarzschild-cells for closed, flat, and open LW 
universes. Each time, a photon was propagated outwards 
in various directions from an initial radius of 10 R, R 
being the Schwarzschild radius. The initial direction of 
travel was given, in terms of block co-ordinates, by 

k a = (1, cos0 n , sin0„), 

with the range of 9 n starting from the purely tangential 
direction of 6q = 7t/ 2 and decreasing until the direction 
was almost completely radial. Both LW and FLRW red- 
shift factors were computed whenever the photon, while 
travelling outwards again, passed an observer co-moving 
with the source. 

We can in principle re-scale the initial k by a con¬ 
stant factor A, but this should not affect the redshifts: 
such a re-scaling would not alter the corresponding null 
geodesic, so the space-time points at which the observer 
and photon meet would remain unchanged; at these 
points, the only change in the photon’s frequencies would 
be a re-scaling by A, but this factor would then cancel out 
of the redshifts. We verified this invariance of the red¬ 
shifts by simulating a photon with different scalings of 
k, and the results, not shown, were indeed invariant for 
at least 11 significant figures; any discrepancies can be 
attributed to numerical error. 

All length-scales in our simulations have been specified 
in terms of R. By simply re-scaling R in one set of results, 
we can readily obtain the results for an equivalent sim¬ 
ulation where the only difference is the magnitude of R. 
In particular, because redshifts are dimensionless quan¬ 
tities, they would not depend on the choice of R, so we 
therefore made the arbitrary choice of setting R = 1 for 
all our simulations. 

We have chosen the dimensions of our Regge block 
as follows. The angular length 2^ was chosen so that 
A 4> = 27t/( 3 x 10 10 ). For the closed universe, the ra¬ 
dial length 2di was chosen to represent a fixed interval of 
Ar = 10 -7 R. For the flat and open universes, it was in¬ 
stead chosen to lengthen for blocks further away from the 
cell centre. This lengthening was implemented so as to 
increase computation speed with only a marginal expense 
to the accuracy: the underlying continuum Schwarzschild 
space-time becomes flatter as one moves further away 
from the centre, so a flat Regge block would approximate 
the region more accurately. Our exact method for deter¬ 
mining the block’s length has been described in Appendix 
ICl Although Ar is no longer constant, it is constrained to 
be an integral multiple of a minimum interval Aro, a pa¬ 
rameter we can freely specify, and the Regge grid is fixed 
so that its first set of blocks covers the region between 
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FIG. 13. A plot of Zlw against Zflrw for 15 trajectories in a flat universe with an initial cell size of rb 0 = 3 x 10 4 R. The 
initial angle 6 n of the trajectory is given in the legend. Each trajectory was traced across 50 cells. A linear regression has 
been performed for each graph, with the first five data points excluded so as to focus only on the linear regime. The regression 
equations and corresponding root-mean-squares of the residuals are listed above in order of decreasing d n . 
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FIG. 14. A plot of redshifts 2 against the radius of observation r 0 b a for photons travelling in the flat universe with an initial 
cell size of rb 0 = 3 x 10 4 R. The photons were travelling initially in the directions of $o = tt/ 2 and 029 = 31l7r/3000. Plots for 
both LW and FLRW redshifts are shown. All four series start at the same data point, and there is a clear jump in the LW 
graphs between the first and second data points. 
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FIG. 15. A plot of redshifts 2 against the radius of observation 
r 0 bs for a photon starting instead with a position and trajec¬ 
tory corresponding to the second data point of the do = 7t/2 
trajectory in Fig. [TJ] Plots for both LW and FLRW redshifts 
are shown. This time, both graphs extend smoothly out from 
the zero-redshift data point. 


R + Aro and R + 2Aro 0 As we shall see below, Aro 
was chosen according to which universe was being sim- 


10 As an informal check on the accuracy, we simulated photons 
propagating across 100 cells in the flat, E = 1 universe us¬ 
ing both fixed and increasing block-lengths. The specific sim¬ 
ulation parameters were an initial Schwarzschild-cell radius of 
rv, (i = 3 X 10 4 R with an initial photon radius of 1 OR: the 


ulated. Finally, the temporal interval was chosen to be 
At = 10 Aro for the flat and open universes and simply 
At = 10 Ar for the closed universe. 

In this section we shall present the redshift results of 
our simulations. We begin with the flat universe, for 
which £1 = 1. To investigate the effect of the initial cell 
size rb 0 on the redshifts, we have simulated flat universes 
for a range of initial sizes from rb 0 = 3 x 10 4 R to rb 0 = 
10 8 R. Following the example of Clifton and Ferreira, we 
have chosen the largest initial size to approximate cells 
with Milky Way-like masses at their centres, as this is 
thought to best represent our actual universe. Depending 
on the initial cell size, Aro ranged from Aro = 10 -5 R 
for rb 0 = 3 x 10 4 R to Ar 0 = 10 -3 R for rb 0 = 10 8 R. 

Figure [13] plots z LW against z F lrw for the universe 
where n, 0 = 3 x 10 4 R. Each trajectory was traced across 
50 cells; only results for 15 trajectories are shown, al¬ 
though we simulated trajectories for 30 angles 9 n rang¬ 
ing from 0o = 7r/2 to 029 = 31l7r/3000 in decrements of 
4l7r/3000. Apart from a brief curve at the start, each 
graph clearly demonstrates a strong linear relationship 
between z LW and z F lrw, and the gradient of the line is 
different for each angle. There are several differences, 
however, between our results and those of Clifton and 
Ferreira. Clifton and Ferreira’s graphs showed more ini- 


Regge blocks’ non-radial dimensions were At = 10 Ar and 
A </> = 27r/(3 x 10 7 ); the radial dimension Ar was fixed at 10 ~^R 
for the fixed-length simulation and was set to be at least that 
length for the increasing-length simulation. The two simulations’ 
redshift results, not shown, agreed for at least five significant fig¬ 
ures, even at the largest radii; however the fixed block-length 
simulation required about 6.5 hours of computation time, while 
the increasing block-length simulation required only 20 minutes. 
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FIG. 16. A plot of photon frequencies v LW against the radius of observation r 0 bs for a photon travelling initially in the 
direction of 0o = tt/ 2; the initial cell size is rb 0 = 3 x 10 4 R. Each graph represents the photon’s frequencies within a single 
cell; these are the frequencies that would be seen by a co-moving observer if the photon intercepted the observer at r 0 bs • 
The analytic frequencies are given by CD, and graphs for the first seven cells traversed are shown. Frequencies were also 
computed numerically by simulating the propagation of a photon across the first cell only; the corresponding graph is shown 
and completely overlaps with its analytic counterpart. 
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tial scatter, which varied depending on the trajectory’s 
initial angle, but their graphs would eventually con¬ 
verge upon a common mean graph as the photon passed 
through an increasing number of cells. Our graphs do not 
display such scatter nor any common mean. However, 
their common mean agreed rather closely with FLRW 
redshifts, as it was given by 1 + z LW ~ (1 + z FLRW )°- 98 
[4ll |: thus in their model, z LW and z FLRW are almost com¬ 
pletely linear. Our graphs, apart from the initial curve, 
are completely linear, but the gradients depend on the 
trajectory followed. 

However Clifton and Ferreira’s relation has the desir¬ 
able feature of passing through the origin, since there 
must be zero redshift at the start of any trajectory; that 
is, both z LW and z FLRW must be zero. All our graphs in¬ 
stead show a jump between the origin and the next data 
point, and this was true of our simulations of all other 
universes as well. When we plotted the redshifts against 
the radius of observation for the rb 0 = 3 x 10 4 R flat 
universe, as shown in Fig. [lTJ we found that the z FLRW 
graph would extend naturally outwards from the zero- 
redshift point at the starting radius of 10 R, but that the 
z LW graph would jump suddenly from the zero-redshift 
point to the next data point. Suppose we ignored our 
current zero-redshift data point and assumed the photon 
actually began its trajectory at the next data point; then 
when we re-calculated all subsequent redshifts based on 
our new initial frequency for the photon, we found that 
the resulting z LW graph progresses naturally from the 
new zero-redshift point to the next data point without 
any jumps, as shown in Fig. H50 

To investigate this initial jump further, we note that 
in universes using CF co-ordinates, the frequency v LW 
measured by a co-moving observer can be expressed as 
a function of the radius at which the measurement is 
made, that is, the radius r 0 b s at which the photon and 
observer intercept. This implies the redshifts z LW can 
be expressed as a function of r 0 b s • Recall that the 4- 
velocity of a co-moving observer is given by (1631) in Regge 
Schwarzschild space-time, with E = Eb and J = 0; this 
vector is a function of the observer’s radial position alone. 
The 4-momentum of the photon is given by the same 
relation as well but with E = E p h and J an arbitrary 
constant; as long as the photon does not cross into the 
next cell, E p h will be constant, so this vector is also a 
function of the photon’s radial position alone. When the 
observer and photon intercept, they will have the same 
radial position, that is, the same value for r, and we 
have denoted above this common r by r 0 b s . By taking 
the scalar product of these two vectors, we obtain the 
measured frequency as a function depending only on the 


value of r 0 bs , 


i/ Eb E ph 

Vlw ~ l 2m 

'f'obs 

E p h _ 

_ 2m j,2 

robs ° ba 

(70) 




We note that if the photon were to cross into the next 
cell, this relation would still hold, but from the bound¬ 
ary conditions, E p h only would change 0 In Fig. [HU and 
Fig. [T71 we have plotted this function for photons travel¬ 
ling at initial angles of n/2 and 31l7r/3000, respectively, 
in the rb 0 = 3 x 10 4 i? flat universe. We have included 
plots for different E p h corresponding to the first few cells 
traversed. We have also simulated a photon’s propaga¬ 
tion across the first cell for both trajectories, and for a 
selection of radii, we have computed the frequencies that 
a co-moving observer would measure if at those radii. 
These numerical results are also included in the figures, 
and we see that they agree very closely with their an¬ 
alytic counterparts. In Fig. [18l we show redshifts z LW 
instead against r 0 b s for the 0q = 7r/2 photon travelling 
through cell 1 only. 

Based on these graphs, we make some speculations on 
the origin of the jump as well as the initial curve in the 
z LW vs z FLRW graphs. We note that for large r Q b s , the 
frequency asymptotically approaches 


= s/E p h (^\Z~Eb — \/Eb — 1 ^ 
= \jE ph for E h = 1. 


Equivalently, the redshift asymptotically approaches 

Zasvm = y/E^{y/m-yfw= ry" 1 

= — 7= ~ 1 for E b = 1. 

V ^ph 

As the photon traversed ever more cells in our simula¬ 
tions, we found that it would intercept the co-moving 
observer at ever larger r 0 b s . So after enough cells, z LW 
should be very close to the asymptotic value z asyrn . Also 
from our simulations, we found that E p h would decrease 
with each subsequent cell implying that z asyrn would in¬ 
crease, which is consistent with the z LW vs z FLRW graphs’ 
positive gradients. Thus the asymptotic, large-radius be- 


11 We did check this result by simulating a photon starting at the 
same radius and direction as that of the original photon when it 
generated the second data point; the resulting graph was indeed 
identical to that of Fig. 1 15 1 


12 Although we have an analytic function for zlw based on r 0 b s , 
we must rely on the simulation to tell us where r 0 & s is, that is, 
where the photon intercepts the observer. And although we can 
also determine E p h analytically from the boundary condition, to 
do this, we still need to know the radius at which the photon 
and cell boundary intercept, and we must rely on the simulation 
to tell us this radius as well. 
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TABLE I. Linear regression equations and the root-mean-squares of the corresponding residuals for Zlw versus Zflrw graphs 
from simulations of photons in the flat universe. The initial cell size was rb 0 = 3 x 10 4 R. This is the same simulation as for 
Fig. 1191 The constants in the regression equations here are much smaller than those in the regression equations in Fig. [13] 


Initial angle 


RMS of residuals 


1.5708 

Zlw 

= 

0.9973 z FLRW 

+ 

4.75 

X 

rc 

1 

o 

t-H 

1.84 

X 

10 -b 

1.4849 

Zlw 

= 

0.9966 Zflrw 

+ 

4.07 

X 

10 -3 

1.88 

X 

10" 5 

1.3991 

Zlw 

= 

0.9959 Zflrw 

+ 

3.41 

X 

10" 3 

2.03 

X 

10 -5 

1.3132 

Zlw 

= 

0.9952 z FLRW 

+ 

2.78 

X 

10" 3 

2.29 

X 

10~ 5 

1.2273 

Zlw 

= 

0.9944 z FLRW 

+ 

2.18 

X 

10" 3 

2.63 

X 

10" 5 

1.1414 

z lw 

= 

0.9937 z FLRW 

+ 

1.62 

X 

10" 3 

3.05 

X 

10~ 5 

1.0556 

Zlw 

= 

0.9929 z FLRW 

+ 

1.09 

X 

10" 3 

3.53 

X 

10~ 5 

0.9697 

Zlw 

= 

0.9922 z FLRW 

+ 

6.06 

X 

10 -4 

4.04 

X 

10“ 5 

0.8838 

Zlw 

= 

0.9915 z FFR w 

+ 

1.75 

X 

10 -4 

4.89 

X 

10 -5 

0.7980 

Zlw 

= 

0.9908 z FLRW 

— 

2.35 

X 

10" 4 

5.49 

X 

10" 5 

0.7121 

Zlw 

= 

0.9902 z FLRW 

— 

6.06 

X 

10 -4 

6.07 

X 

10~ 5 

0.6262 

Zlw 

= 

0.9897 z FLRW 

— 

9.37 

X 

10 -4 

6.62 

X 

10~ 5 

0.5404 

Zlw 

= 

0.9892 z FLRW 

— 

1.23 

X 

10" 3 

7.12 

X 

10~ 5 

0.4545 

Zlw 

= 

0.9888 z FLRW 

— 

1.48 

X 

10 -3 

7.57 

X 

10~ 5 

0.3686 

Zlw 

= 

0.9884 z FLRW 

- 

1.69 

X 

10" 3 

7.94 

X 

10” 5 


TABLE II. Linear regression equations and the root-mean-squares of the corresponding residuals for Zlw versus z F lrw graphs 
from simulations of photons in the rb 0 = 10 8 R flat universe. This is the same simulation as for Fig. l21a| and Fig. l21bl Regression 
was performed on the z F lrw > 0.001 region. Both the regression equation constants and the residuals are even smaller here 
than in Table jT] 


Initial angle 

L5708 

1.4635 

1.3561 

1.2488 

1.1414 

1.0341 

0.9268 

0.8194 

0.7121 

0.6048 

0.4974 

0.3901 


Regression equation 

z LW = 1.0215 z FLRW + 5.4 x 10~ b 
Zlw = 1.0200 Z FLRW + 2.9 X 10 -5 
= 1.0167 Z FLRW + 1-5 X 10“ 5 
Zlw = 1.0130 Zflrw + 8.8 X 10~ 6 
= 1.0095 z FLRW + 5.4 x 10" 6 
= 1.0066 Z FLRW + 3.4 X 10" 6 
= 1.0041 z FLRW + 2.2 x 10" 6 
= 1.0020 Z FLRW + 1-4 X 10~ 6 
Zlw = 1*0002 Z F lrw + 9.2 X 10~ 7 
Zlw = 0.9988 Zflrw + 4.5 X 10 - ' 
Z ivv = 0.9976 Zflrw + 1-0 X 10" 7 
Zlw = 0.9967 Zflrw - 1-6 X 10~ 7 


RMS of residuals 

1.66 x 10~ (> 

1.45 x 10" 6 
1.04 x 10" 6 
6.97 x 10 -7 
4.62 x 10- 7 
3.08 x 10~ 7 
2.04 x 10 -7 
1.31 x 10" 7 
6.34 x 10~ 8 
3.10 x 10~ 8 
7.08 x 10~ 9 
1.24 x 10 -8 


haviour of the redshifts is responsible for the linear be¬ 
haviour of the z LW vs z FLR w graphs. The graphs’ ini¬ 
tial behaviour however must be explained by the be¬ 
haviour of the redshifts at low r a b Sl where the photon 
is still sufficiently close to the central mass to feel its 
influence strongly. We note that the frequencies always 
approach Vasym from below, implying that the redshifts 
always approach z asyrn from above. This is consistent 
with the initial curve in the z LW vs z FLRW graphs, as this 
curve converges into the linear regime from above as well. 
At much smaller r 0 b s however, the redshift may deviate 
much more significantly from the asymptotic value. In 
particular, the deviation of the zero-redshift point, where 
r 0 b s is smallest, may be much greater than the deviation 
of the next data point, where r 0 b s may have increased suf¬ 
ficiently for the redshift to be much closer to the asymp¬ 
totic value. This would account for the sudden initial 


jump seen in the z L w vs z FL rw graphs. Therefore, the 
initial jump and the initial curve in the z L w vs z F lrw 
graphs can be attributed to the photon being closer to 
the central mass and thus feeling its effects more strongly. 

To test our conjecture, we simulated photons in the 
same universe but starting from a large radius, specif¬ 
ically r = 10 4 R, such that the central mass’ influence 
would be much reduced. Fig. ll9al shows the resulting plot 
of z L w against z FLRW for the most radial and the most 
tangential trajectories simulated. The graphs now follow 
completely straight lines passing very close to the ori¬ 
gin, with the initial jump completely gone and the initial 
curve nearly absent. All other trajectories at intermedi¬ 
ate angles displayed similar behaviour, with the graphs’ 
gradients increasing as the initial trajectory became more 
tangential. We also note that the gradients are nearly 
unity, indicating that z L w is nearly identical to z F lrw■ 














22 



+ cell 1 (E p h = 0.9), numerical 
X cell 1 (E p h = 0.9), analytic 
• cell 2 (E p h = 0.8795), analytic 
▼ cell 3 (E p h = 0.8596), analytic 


r obs/ R 

▲ cell 4 (Eph = 0.8403), analytic 

♦ cell 5 (Eph = 0.8214), analytic 

# cell 6 (Eph = 0.8032), analytic 
A cell 7 (Eph = 0.7854), analytic 


FIG. 17. The equivalent plot to Fig. 1161 but for a photon travelling initially in the direction of O 29 = 31l7r/3000. Once again, 
the numerical and analytic graphs for cell 1 overlap completely. 
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r obs/R 


+ cell 1 (Eph = 0.9), numerical 


X cell 1 (Eph = 0.9), analytic 


FIG. 18. The equivalent plot to Fig. [16] but showing the redshifts Zlw of cell 1 only against the radius of observation r 0 bs • The 
two graphs overlap each other completely. 


In Fig. I19bl we have plotted the difference between the 
z LW graphs and the reference graph of z LW = z FLRW . 
The deviation from linearity at low r Q b s becomes clearer 
now, although we also see that it is short-lived. From 
this figure, we see that linearity is relatively well estab¬ 
lished by z FLRW w 0.2; thus for each trajectory, we have 
regressed the corresponding z LW vs z FLRW graph for the 
region z FLRW > 0.2. All regression equations and cor¬ 
responding root-mean-squares of the residuals have been 
listed in Table [TJ and two of these regressions are shown 
in Fig. I19al These equations show that our observations 
for the graphs in Fig, llflal holds generally for all trajecto¬ 
ries, specifically that the gradients are nearly unity and 
that the graphs pass very close to the origin. Thus, the 
influence of the masses on redshifts has been significantly 
suppressed, and the redshifts now behave nearly identi¬ 
cally to FLRW, much like Clifton and Ferreira’s mean 
redshift result. 

For flat universes with larger initial cell sizes, our sim¬ 
ulations yielded results similar to those of the rf, 0 = 
3 x 10 4 R universe. Figure [20] depicts plots of z LW against 
z FLRW for the do = 7t/2 trajectory for universes with rb 0 
ranging from 10 5 to 10 8 . All graphs again converged 
quickly to a straight line from above, including graphs, 
not shown, for the other trajectories. One noticeable 
development though was that the curve at the start of 
the graphs became more pronounced as the initial size 
increased. However when we simulated photons starting 
from a large radius again, the initial curve was again sup¬ 
pressed, and the initial jump was again absent, as shown 
in Fig. I21al Fig. 12lbt and Table [TT1 for photons starting 
from r = 10' I? in the n, 0 = 10 s A? universe; the initial 
curve is even more short-lived, as seen in Fig. 12 lbl and 
the regression lines for all trajectories pass even closer to 


the origin, as demonstrated by the even smaller constants 
in the regression equations in Table [II] 

We also notice, in Fig. [131 that the gradients of z LW 
vs z FLRW graphs decrease with 9 n . In fact, if we plot the 
gradients against cos(0„), we find a strongly linear rela¬ 
tionship between the two quantities, as shown in Fig. 1221 
As the figure shows, this is also true of the gradients for 
all the other flat universes we simulated. Since cos (9 n ) 
determines the radial component pi n it of the photon’s 
initial velocity, this implies that the gradients are lin¬ 
early related to pinit /fi n it . As the photon’s initial veloc¬ 
ity becomes increasingly radial, the photon’s frequency 
becomes increasingly blueshifted relative to the FLRW 
redshift. We believe this is due to the central mass’ 
stronger influence on the photon, as a more radial trajec¬ 
tory would bring the photon closer to the central mass. 
Quite remarkably, this influence can alter the lattice uni¬ 
verse redshifts z LW from their FLRW counterparts z FLRW 
by as much as 30%, as shown by the right end of the 
rb 0 = 10 8 graph. Since the lattice universe can generate 
lower redshift values, an observer fitting these redshifts 
to an FLRW model could significantly underestimate the 
age of the photons’ source. 

Everything we noticed about redshifts in the flat uni¬ 
verse also applied to open universes. We simulated open 
universes for values of Eb in the range of 1.1 > Eb > 1. 
In each simulation, the initial cell size was fixed to be 
rb 0 = 3 x 10 4 , and photons were propagated along 12 
trajectories with initial angles ranging from 9q = 7r/2 to 
0ii = 1497 t/ 1200 in decrements of 4l7r/1200. Depending 
on Eb, Aro ranged from 10 -5 R for Eb = 1 to 10 -3 R 
for Eb = 1.1. Figure [53] shows a clear linear relationship 
again between z LW and z FLRW for the Eb = 1.1 universe. 
The regression equations’ non-zero intercepts indicate the 
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(a) 



ZFLRW 


1.5708 -t- 

1.4849 -*■ 

1.3991 - *■ 

1.3132 -^ 


1.2273 


0.8838 - 

— 0 - 

0.5404 

1.1414 


0.7980 - 

— V — 

0.4545 

1.0556 


0.7121 

— V — 

0.3686 

0.9697 

—A- 

0.6262 

— V — 



(b) 

FIG. 19. (a) A plot of z LW against zflrw for a tangential and a nearly radial trajectory in the flat universe with initial cell size 
of ri 0 = 3 x 10 4 R. This is the equivalent plot to that of Fig. [lil] but with the photon starting at 10 4 R instead. The regression 
lines now pass very close to the origin, and the initial curve is nearly absent; we also note that the gradients are nearly unity. 
Other trajectories, not shown, showed similar behaviour, (b) A plot of Zlw — Zflrw against Zflrw for 15 trajectories in the 
same universe, including the two trajectories in (a). Each graph still shows an initial curve, but it is small and short-lived, and 
linearity is soon established. 
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Zflrw 


Zflrw 


(a) 


(b) 




Zflrw 


Zflrw 


(c) 


(d) 


FIG. 20. A plot of Zlw against Zflrw for the do = 7r/2 photon trajectory in the (a) rb 0 = 10 5 R universe, (b) r& 0 — 10 6 R 
universe, (c) r(, 0 = 10 7 R universe, and (d) rb 0 = 10 8 R universe. Aro was 10~ 5 R for (a) and (b), 10~ 4 R for (c), and 10~ 3 R 
for (d). 
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Zflrw 

1.5708 + 0.3901 x 

(a) 



Zflrw 


1.5708 -f- 

1.4635 - *■ 

1.3561 - •- 


1.2488 

1.1414 

1.0341 


0.9268 - +- 

0.8194 — At 
0.7121 ©- 


0.6048 — v 
0.4974 — v 
0.3901 — v- 


(b) 

FIG. 21. (a) A plot of Zlw against Zflrw for a tangential and a nearly radial trajectory starting at radius 10' R in the 

Tb 0 = 10 8 R flat universe. Each trajectory was traced across 15 cells. Compared to Fig. I19al the regression lines here pass even 
closer to the origin; the initial curve is even more suppressed; and the gradients are even closer to unity. Other trajectories, not 
shown, showed similar behaviour, (b) A plot of z LW — z F lrw against z F lrw for 15 trajectories in the same universe, including 
the two trajectories in (a). Compared to Fig. I19bl the initial curve here is even smaller and more short-lived. 























gradient of Zlw vs zflrw 


1.20 


1.10 


1.00 - 


0.90 - 


gradient = 
RMS of 
gradient = 
RMS of 
gradient = 
RMS of 
gradient = 
RMS of 
gradient = 
RMS of 


= -0.3477 cos(0) + 
residuals = 7.01 X 
= -0.3380 cos(0) + 
residuals = 2.11 X 
= -0.3322 cos(0) 
residuals = 2.16 x 
= —0.3234 cos(0) 
residuals = 2.11 X 
= -0.3104 cos(0) + 
residuals = 1.99 X 




0.80 - 

+ r b = 3 X 10 4 
X r b = 10 s 
0.70 |- • r b = 10' 

T n = 10 7 
A r b = 10 8 

0.60 


0.00 


0.10 


0.20 


0.30 
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0.90 


1.00 


FIG. 22. A plot of the gradients of Zlw vs Zflrw graphs against cos^n) for flat universes of different initial cell size r bo . 
linear regression has been performed on each graph, and the regression equations are displayed in order of increasing r bfJ . 
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presence again of a jump between the zero-redshift point 
and the first non-zero redshift point. Similar behaviour 
was seen for universes corresponding to other values of 
£/,. Figure [23] shows that for all Eb, there is again a clear 
negative linear relation between the gradients of z LW vs 
Zflrw graphs and cos (0 n ), or equivalently between the 
gradients and Pinit/nnit- 

To see how the gradients of z LW vs z FLRW change with 
Eb, we have plotted the gradients against Eb in Fig. [25] 
for a selection of trajectories. The graphs indicate that 
the gradients approach some asymptotic value as Eb in¬ 
creases and that this asymptotic gradient decreases as 
the photon trajectory becomes more radial, consistent 
with the behaviour illustrated in Fig. [2~il 

Finally, we shall now discuss our simulation of the 
closed universe. We have chosen to focus on an LW uni¬ 
verse built from the 600-cell Coxeter lattice described in 
Appendix E as it is the most finely subdivided closed 
Coxeter lattice possible. Using (1551) and (1351) . we find 
that Eb ~ 0.96080152145 for this universe, and the max¬ 
imum cell size is therefore approximately 25.51?. If we 
embed a Schwarzschild-cell into the comparison hyper¬ 
sphere, as described in Section m then the ratio be¬ 
tween the angular distance, il>h, from the cell centre to 
the horizon, as measured from the hypersphere’s centre, 
and the corresponding distance, ipb, to the cell bound¬ 
ary is approximately 0.0392(3 for comparison, Clifton et 
al. [13 have shown that on the time-symmetric hyper¬ 
surface, the exact 3-metric is conformally equivalent to 
that of a 3-sphere, and on this 3-sphere, the equivalent 
ratio is 0.0147; our value is about 2.6 times greater than 
theirs, although their ratio uses the shortest angular dis¬ 
tance to the cell boundary, which is not spherical in the 
exact lattice. In any case, the 600-cell universe is much 
smaller than the universes we have heretofore been con¬ 
sidering, and the central mass’ influence will therefore be 
much stronger. We note that in the LW co-ordinates for 
the closed universe, we do not have freedom to choose 
the initial size of the cell; once the photon’s initial co¬ 
ordinates are chosen, then the cell boundary’s position 
is determined, since both the boundary and the source 
must reach their respective maximum radii at the same 
time t. For our simulation of this universe, we chose 
to propagate photons along 30 trajectories ranging from 
00 = to 029 — going in decrements of 

The simulation results were again similar to those for 
the flat and open universes but with several slight differ¬ 
ences this time. The graphs of z LW against z FLRW again 
followed straight lines but the gradient was slightly dif¬ 
ferent for when the photon was outgoing and when it was 
ingoing, as shown in Fig. [26] As the figure also shows, 
the graphs did not necessarily intercept the origin either. 
Several of the graphs for ingoing photons also had a sub- 


13 This value has been calculated using dD , 1371 ). 1231 ). and (1391 ; 
from these equations, it follows that sin 0 ^ — sin 3 - 0 b, where 0 b 
is given by gsj for TV = 600. 


tie bend in their low z end. We believe this bend to be 
analogous to the initial curve we saw in the flat universe 
graphs previously: it corresponds to redshifts measured 
at low r 0 b s and is more pronounced in the more radial tra¬ 
jectories, that is, the smaller 9 n trajectories, which pass 
closer to the central mass; thus, like the initial curve in 
the flat universe graphs, we believe the bend is caused by 
the central mass’ stronger influence at low r 0 b s (3 In Fig. 
[271 we show graphs of z LW against z FLRW for when the 
photons are outgoing, and in Fig. [28l we show the corre¬ 
sponding graphs for when the same photons are ingoing. 

In Fig. [251 we have plotted the gradients of z LW vs 
Zflrw against cos(0„). Two graphs are shown, one for 
the gradients when the photons are outgoing, and the 
other for when they are ingoing. To exclude the initial 
bend in the graphs for ingoing photons, we excluded data 
points corresponding to r 0 b s less than 10 R, the photons’ 
starting radius(3 The figure shows that, as with the 
other universes, the gradient decreases as the photon’s 
trajectory becomes more radial, but now, the relationship 
between the gradients and cos(0„) is no longer linear. 

Thus to summarise, our simulations have demon¬ 
strated several features common to redshifts in all LW 
universes. The LW redshifts generally increase linearly 
with their FLRW counterparts in any universe and for 
any photon trajectory. When r 0 b s is small such that the 
central mass has a stronger influence, there is some devia¬ 
tion away from this linear behaviour. Any influence from 
the central mass can be suppressed by starting the photon 
at very large radii; z LW then becomes completely propor¬ 
tional to z FLRW . The LW redshifts also generally decrease 
relative to their FLRW counterparts as the photon takes 
a more radial trajectory; this can also be attributed to 
the stronger influence of the central mass since a more 
radial trajectory would pass closer to it. The mass’ in¬ 
fluence on redshifts can sometimes be quite significant, 
causing LW redshifts to deviate from their FLRW coun¬ 
terparts by as much as 30%. This would have significant 


14 We note that the earlier analysis of the dependence of ulw on 
r obs> given by function 1701) . would require some modification in 
this context. In particular, as we are now using LW co-ordinates, 
function 1701) would change, with parameter E 5 replaced by a 
parameter E 0 corresponding to the constant E for the co-moving 
observer’s geodesic. The new function though would still have 
the same functional form as the original and therefore the same 
general behaviour. However since the radius of the observer’s 
geodesic is bounded, the new function’s domain is bounded from 
above by r 0 b s < 2m/(1 — E 0 ), otherwise the square-root in the 
function would turn imaginary; therefore there cannot be any 
large-r 0 ;, s asymptotic regime in the new function, and hence this 
aspect of our previous analysis is not transferable to the closed 
universe. Nevertheless, we believe that the general conclusion 
still applies, that the initial jump and the bend in zlw vs zflrw 
graphs is caused by the central mass’ stronger influence at low 

^obs • 

15 In fact, including all data points in the regression introduced 
some ‘jaggedness’ in the graph of ingoing gradients, especially 
towards the low cos(0 n ) regime. By excluding these points, the 
‘jaggedness’ smoothed away, and we were left with a curve much 
more similar to that for the outgoing gradients. 
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FIG. 23. A plot of Zlw against z F lrw for 12 trajectories in the Eb = 1.1 universe. For this simulation, Aro was 1CP 3 R. The 
initial angle 9 n of the trajectory is given in the legend. Each trajectory was traced across 15 cells. A linear regression has been 
performed for each graph. The regression equations and corresponding root-mean-squares of the residuals are listed above in 
order of decreasing 9 n . 
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FIG. 24. Plot of the gradients of Zlw vs zflrw graphs against cos (9 n ) for open universes of different values of E b . Also shown 
is a plot for the flat universe, corresponding to E b = 1. A linear regression has been performed on each graph, and both the 
regression equation and the root-mean-square of the residuals are displayed in order of increasing E b . 
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FIG. 25. A plot of the z LW vs z FLRW gradients against E b 
for a selection of initial photon trajectories. Each graph ap¬ 
proaches an asymptotic value as E b increases. 

implications if we attempted to fit redshift data from a 
lattice universe onto an FLRW model, as quantities such 
as the age of the universe may be incorrectly estimated. 



Zflrw 

FIG. 26. A plot of Zlw against Zflrw for two trajectories 
of photons, 9q — 7r/2 and 029 = 177 t/ 150. Each trajectory’s 
graph is divided into a subgraph for when the trajectory is 
outgoing and a separate subgraph for when the same trajec¬ 
tory becomes ingoing; the two subgraphs meet at the redshift 
data point closest to maximum expansion of the universe. 
There is a slight bend in the blueshift regime of the ingoing 
029 = 177r/150 graph. 

Thus we see certain significant effects arising from the 
dumpiness’ of the LW universe which an FLRW-based 
model could not adequately capture. 
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FIG. 27. A plot of Zlw against Zflrw while photons are outgoing. Graphs for a selection of photon trajectories 9 n are shown. 
A linear regression was performed on each graph, and the regression equation is listed in order of decreasing 9„. 



32 


This influence of the masses on redshifts can be un¬ 
derstood as an integrated Sachs-Wolfe effect (ISW). The 
ISW describes the net redshift induced on a photon as 
it passes through a fluctuating gravitational potential 
caused by fluctuations in the energy density [4814511] . In 
a flat FLRW universe, a time-varying potential is nec¬ 
essary to generate a non-zero ISW; however when the 
universe is matter-dominated, the potential is static, and 
therefore the ISW would be zero. If there is a non-zero 
cosmological constant though, the potential in a matter- 
dominated universe does become time-varying, and hence 
a non-zero ISW can result [ 52 }. Indeed, recent precision 
observations of the cosmos have not only established that 
our universe is very nearly flat but that there is definitely 
a non-zero ISW [&$, HH ; naturally, this has been inter¬ 
preted as evidence of a non-zero cosmological constant 
in our universe. However, our results suggest that this 
might be explainable within a matter-dominated universe 
without needing a cosmological constant; we may simply 
require a model that better reflects the inhomogeneous 
matter distribution than FLRW space-times allow. Be¬ 
cause of the significant implications a non-zero ISW may 
have on the dark energy question, this modelling problem 
deserves further investigation. 

VIII. DISCUSSION 

We have investigated the properties of the Lindquist- 
Wheeler universes, as we hope they might provide some 
insight into what observable effects the ‘lumpy’ matter 
distribution of the actual universe might yield. And al¬ 
though the LW universes are only approximations rather 
than exact solutions to the Einstein equations, we be¬ 
lieve they model enough of the underlying physics to 
yield at least meaningful qualitative insights into the be¬ 
haviour of the actual universe. Much of the LW uni¬ 
verses’ dynamics bear strong resemblance to those of the 
matter-dominated FLRW universes. Additionally, pho¬ 
ton redshifts in LW universes behaved roughly similarly 
to their FLRW counterparts. Yet there were also cer¬ 
tain direction-dependent effects in the redshifts due to 
the ‘lumpiness’ of the universe’s matter distribution. 

We noted that there were some differences between 
Clifton and Ferreira’s results and ours for the flat uni¬ 
verse. Most notably, our redshifts displayed the effects 
of the inhomogeneous matter distribution more. The ef¬ 
fects became stronger the closer the photon passed to a 
mass and the more radial the photon’s trajectory was. 
As a result, our redshifts demonstrated a non-zero ISW, 
which was not present in Clifton and Ferreira’s results. 
Otherwise, when the photon was far from any mass, its 
redshift would behave much more similarly to FLRW 
redshifts as well as to Clifton and Ferreira’s mean red¬ 
shift relation. The main difference between Clifton and 
Ferreira’s model and ours was the definition of the cell 
boundary and the choice of conditions for propagating 
particles across boundaries. As Clifton and Ferreira have 


shown, these can significantly alter the redshift relations; 
for instance, using a different choice of boundary condi¬ 
tions, Clifton and Ferreira obtained the redshift relation 
of U -z LW « (1 +z FLRW ) 7 / w instead [13,1111, which clearly 
agreed far less with FLRW redshifts. Further investiga¬ 
tion is required to determine whether our discrepancy 
from Clifton and Ferreira’s results is indeed significant 
or merely an artefact of our boundary conditions which 
somehow exaggerates the effects of the masses. Neverthe¬ 
less, neither our inhomogeneous results nor Clifton and 
Ferreira’s more homogeneous one can be ruled out at the 
moment. 

There are several ways in which our investigation can 
be extended. It would be interesting to examine the op¬ 
tical properties of the LW universe, as this may have 
important implications for actual astrophysical observa¬ 
tions. Clifton and Ferreira have already done this for 
their implementation of the LW universe, but we have 
used a different implementation from theirs, and this has 
led to certain differences in the behaviour of redshifts. 
However, Clifton and Ferreira have suggested that the op¬ 
tical properties, unlike the redshifts, may actually be in¬ 
sensitive to the choice of boundary conditions [4J|. Nev¬ 
ertheless, it would be very interesting to see whether that 
is indeed the case for our boundary conditions. It would 
also be interesting to extend our study to LW universes 
with a non-zero cosmological constant and attempt to 
evaluate by how much inhomogeneities reduce the need 
for such a constant. Clifton and Ferreira have already 
constructed an appropriate extension of the LW universe 
based on the Schwarzschild de-Sitter metric, and they 
have shown that the corresponding Friedmann-like equa¬ 
tion strongly resembles its FLRW counterpart as well. 
Moreover, using their boundary conditions, they have 
also shown that the cosmological constant density ratio 
Ua can be reduced by about 10 % [ 4 j], [55J , and again this 
was rather insensitive to which boundary condition they 
used. It should be possible to include Clifton and Fer¬ 
reira’s A-Schwarzschild-cells in our implementation and 
to investigate the resulting model. Finally, our model has 
a still very idealised distribution of matter. Each mass 
is identical and distributed on a perfect lattice which is 
clearly not the case in the actual universe. We should like 
to extend our model to allow for different sized masses 
and cells. To this end, we have derived in Appendix [D] a 
set of conditions that must be satisfied at the boundary 
between two neighbouring cells of different mass m. We 
leave the numerical simulation of such universes and the 
detailed investigation of their properties to future work. 
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FIG. 28. A plot of Zlw against z F lrw while photons are ingoing. Graphs for the same trajectories as in Fig. [23 are shown. 
A linear regression was performed on each graph, and the regression equation is listed in order of decreasing 9 n . Though in 
this case, the regression only includes redshifts measured at r 0 b s equal to at least 10 R, the photon’s starting radius, so as to 
exclude points in the initial bend. 
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TABLE III. All possible lattices obtained by tessellating 3- 
spaces of constant curvature with a single regular polyhedron. 


Elementary 
cell shape 

Number of 
cells at a 
lattice edge 

Background Total cells 
curvature in lattice 

tetrahedron 

3 

+ 

5 

cube 

3 

+ 

8 

tetrahedron 

4 

+ 

16 

octahedron 

3 

+ 

24 

dodecahedron 

3 

+ 

120 

tetrahedron 

5 

+ 

600 

cube 

4 

0 

oo 

cube 

5 

- 

oo 

dodecahedron 

4 

- 

oo 

dodecahedron 

5 

- 

oo 

icosahedron 

3 

- 

oo 


FIG. 29. A plot of zlw vs zflrw gradients against cos(9„), 
where 6 n is the initial direction of the photon. The gradients 
for when the photons are outgoing and for when they are 
ingoing are plotted separately, and all gradients are obtained 
by performing a linear regression on z LW vs Zflrw data. To 
avoid the initial bend in the Zlw vs z F lrw graphs for ingoing 
photons, the ingoing gradients are obtained by regressing over 
only data points for which r Q (, s is at least equal to the photon’s 
starting radius of 10 R. 
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APPENDIX A: REGULAR LATTICES IN 
3-SPACES OF CONSTANT CURVATURE 

In this appendix, we shall list all possible lattices that 
cover 3-spaces of constant curvature with a single regu¬ 
lar polyhedral cell. The cell is tiled to completely cover 
the 3-space without any gaps or overlaps. This tes¬ 
sellation problem has been thoroughly studied by Cox- 
eter (56). Clifton and Ferreira j40] have succinctly sum¬ 
marised Coxeter’s results relevant to our discussion, and 
we have presented their summary in Table m 

Column 1 gives the lattice cell’s shape. Column 2 effec¬ 
tively determines the lattice’s structure: it indicates how 
many of Column l’s elementary cells meet at any lattice 
edge. Column 3 indicates whether the background 3- 


space has positive, flat, or negative curvature. Column 4 
gives the number of cells needed to cover the 3-space; only 
lattices on positively curved space can have a finite num¬ 
ber of cells. To construct a closed lattice universe, one 
has six choices of lattices; for example, one can choose a 
600-cell lattice where cells are equilateral tetrahedra and 
where five tetrahedra meet at any edge. To construct a 
flat lattice, one has but a single choice. To construct an 
open lattice, one has four. 

APPENDIX B: RADIAL VELOCITIES IN 
REGGE SPACE-TIME 

In this appendix, we shall present our numerical results 
supporting (1551) to be the 4-velocity tangent to radial 
time-like geodesics in Regge Schwarzschild space-time. 
These geodesics include, most importantly, those of test 
particles co-moving with a Schwarzschild-cell boundary. 
We began by determining numerically the escape velocity 
for a test particle. We propagated a test particle radially 
outwards from a series of initial radii 7+^ and with a 
series of initial velocities Pi n it- We started pi n n from 0 
and increased it until the maximum radius r m ax attained 
by the particle was very large. For each we then 

plotted R/r max against its corresponding pi n u, where R 
is the Schwarzschild radius. Both rj„j t and r max refer to 
the lower Schwarzschild label of the block in which the 
particle is found. Our plots are shown in Fig. [301 Each 
graph is very well-fitted by a quadratic curve of the form 

R . .2 „ 

- - — - A Pi n it + B. 

<max 
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“I" T init/R — 8 X Tinit /R — 32 # Tin.it/— 128 ▼ T init /— 512 

FIG. 30. A plot of R/rmax versus pinit for various rtnit . Each graph has 100 data points. A quadratic regression has been 
performed on each graph, and the regression equations are ordered from ri n u = 8 R to 512 R. The regressions were performed 
without a linear term; regression was also attempted with a linear term present, but it was found to be many orders of 
magnitude smaller than the other two terms. The simulation’s block parameters were R = 31, At = lOAr, Ar = A/10 5 , and 
A (t> = 2tt/( 3 x 10 7 ). 
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FIG. 31. A plot of {u p f against p 2 for various r ma x < 0, where u p is given by (1581) . Each plot consists of velocities computed 
at 100 different radii along the course of the test particle’s trajectory. A linear regression was performed on each graph, and 
the corresponding equations are ordered from R/rmax = —1/11 to —1/91. The graphs of the regressions completely overlap 
each other. The gradients only begin deviating from unity at the 10 -5 order of magnitude. The simulation’s block parameters 
were Tin.it =2 R, R = 21, At = Ar = R/ 10 s , and Arj> = 2n/(3 x 10 7 ). 


Moreover, the co-efficients A and B always satisfy the 
relations A + B = 1 and B = I?/ r imiEl From these rela¬ 
tions, we can therefore infer a relation for p corresponding 
to that given in (1581) . and the f component follows from 
normalisation. 

We have just provided numerical support for (1551) for 


16 We have also looked at graphs for different R and r; n ;t, not 
shown, and they also conform to this pattern. 


test particles following geodesics where r max > 0, corre¬ 
sponding to closed orbits or orbits at the escape velocity. 
At this point, we conjectured that for open orbits, (1551) 
would still apply but with r max < 0. To test this, we 
propagated a test particle outwards with an initial ve¬ 
locity given by (1551) but for r max < 0, and we examined 
how the velocity evolved with r. Let us denote the ve¬ 
locity of our simulated particle by p. Along the particle’s 
trajectory, we compared p 2 against ( u p ) 2 for the same ra¬ 
dius, where u p is given by the radial component of (1581) . 
Our comparison is presented graphically in Fig. [31] In all 
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FIG. 32. A plot of 1/p versus r/R comparing simulated par¬ 
ticle velocities with u p . For the simulated particles’ graphs 
p = p, while for the u p graph, p = u p . The simulation was run 
for various sizes Ar of the Regge block. For all simulations, 
the particle was propagated from r = 2 R to 3 x 10 4 R. The 
simulation’s other block parameters were R = 1, At = Ar, 
and A <j> = 27t/(3 x 10 7 ). 

cases, the graphs’ gradients were effectively unity and the 
constants effectively zero, thus indicating that the par¬ 
ticle’s velocity does indeed obey (1551) . In Table Q3D we 
provide a partial list of the values of p 2 and ( u p ) 2 as well 
as their percentage difference. The percentage difference 
is very small in all rows, although there appears to be an 
increasing trend with increasing r. We find that our con¬ 
jecture is supported by numerical results, and hence (1581) 
appears valid for all types of radial time-like geodesics FI 
We have just found that the initial velocity pinit of a 
simulated particle must equal at least that of (1551) in 
order for the particle to ‘escape’. We next examined 
the long-term behaviour of the particle’s velocity p as 
the particle propagated outwards from an initial velocity 
equalling at least the escape velocity. We found that p 
would not stay equal to u p as given by (1551) but would 
instead decrease at a slower rate as a function of the par¬ 
ticle’s radius. For example if we started a particle at 
the escape velocity at = 2 R and propagated it out¬ 
wards, then by the time it reached r = 3 x 10 4 R, there 
was a very significant discrepancy between p and u p , as 
shown in Fig. [351 with p being consistently larger than 
u p . This discrepancy reduced if we improved the radial 


17 Again, we have also simulated the cases of r m ax = 
— R , —21 R, —41 R, —61 R , and —81 R, not shown here, and they 
all display the same behaviour as in Fig. 1311 The percentage 
difference seen between p 2 and (u p ) 2 is also small but slowly 
increasing in all r marc that we looked at. 


resolution Ar of our Regge blocks, but it was still present 
even for resolutions as high as Ar = R/ 10 5 . If instead, 
using the same method as previously, we numerically de¬ 
termined the escape velocity for a particle starting at 
finit = 3x 10 4 R, we again obtained the same quadratic 
behaviour between R/r max and p as that in Fig.[30l thus 
indicating the escape velocity should still correspond to 
(1551) . Indeed, this result was obtained even when the ra¬ 
dial resolution was as low as Ar = i?/10, as shown in Fig. 
1331 even though there is a huge discrepancy at this reso¬ 
lution between p and u p in Fig. 1321 Therefore even if a 
particle began at the correct escape velocity for its initial 
radius, it would gradually travel faster than the correct 
escape velocity for its subsequent radii as it propagated 
along its trajectory. 



Pinit 


FIG. 33. A plot of R/rmax versus pinit for n„u = 3x 10 4 R. A 
quadratic regression has been performed and the correspond¬ 
ing equation shown. The simulation’s block parameters were 
R = 1, At = lOAr, Ar = R/ 10, and Arf> = 2tt/(3 x 10 7 ). 

The discrepancy between p and u p can be understood 
by considering the manner by which particles are de¬ 
flected in the Regge model. Suppose a particle enters 
a Regge block of radius r,; by the —pi face and travels at 
the escape velocity as given by (1551) . If the particle then 
enters the next block by the —Pi+i face as well, it will 
have undergone no deflection and hence continue having 
the same velocity as at r*. However, (1551) indicates that 
the escape velocity at ri + i should be smaller than that 
at rj, and hence we see the beginning of a discrepancy. If 
the particle passes through a sizeable number of p faces 
in succession, it would accumulate a significant discrep¬ 
ancy such that the particle’s velocity would be notice¬ 
ably larger than the escape velocity for its radius. This 
explains why the graphs for simulated particles in Fig. 1321 
are always smaller than the graph for u p . Only when the 
particle crosses a r surface does it undergo a change in 
velocity, but this change may not be enough to overcome 
the discrepancy already accumulated. By a similar argu¬ 
ment, we expect the radial velocity of a radially in-going 
particle to decrease less quickly than u p in (1551) . because 
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TABLE IV. A list of p 2 and ( u p ) 2 at various radii, and the percentage difference between them. This data is for i?/r maa , = —91. 


r/R 

P 2 

(■ u p ) 2 

% difference 

2 

1.0219794649 

1.02197802197802 

0.00014118894047312 

102 

0.021001226724 

0.0209988015595972 

0.0115477273528382 

202 

0.016021230625 

0.0160188121184031 

0.0150956356191602 

302 

0.014350362849 

0.0143477802197414 

0.0179969613713785 

402 

0.013512667536 

0.0135101783683129 

0.0184209940820187 

502 

0.013009455481 

0.0130069512260618 

0.0192494985037111 

602 

0.012673806084 

0.0126711894819952 

0.0206457475160902 

702 

0.012433811049 

0.0124312192608481 

0.0208446802168433 

802 

0.012253604416 

0.0122511667237755 

0.0198936748870944 

9202 

0.011101572496 

0.0110988891274903 

0.0241710668530401 

9302 

0.011100308164 

0.0110977080370158 

0.0234239171181774 

9402 

0.011099254609 

0.0110965512181986 

0.02435650768163 

9502 

0.011097990409 

0.0110954196820738 

0.023163895728025 

9602 

0.011096936964 

0.0110943114060724 

0.0236602040377857 

9702 

0.011095883569 

0.0110932263900415 

0.0239474300718204 

9802 

0.011094830224 

0.0110921626114406 

0.0240437438476478 

9902 

0.011093776929 

0.011091121081329 

0.0239399772323361 



FIG. 34. A plot of p versus r/R comparing a simulated 
particle velocity with u p . The simulation began at ri„it = 
3 x 10 4 R, but only data points from r = 2 R to 100 R are 
shown. The simulated velocity is always greater than or 
equal to u p . The simulation’s block parameters were R = 1, 
At — Ar = R , and A <j> = 2n/(3 x 10' ). 

if the particle crosses a +p face with the correct escape 
velocity, it will cross the neighbouring +p face with the 
same unchanged velocity, which would be less negative 
than u p for the same radius; and indeed, this is what we 
see in Fig. [Ml By improving the resolution of our Regge 


model, we slow the growth rate of the discrepancy, as we 
saw in Fig. [Ml because our model becomes a better ap¬ 
proximation to the underlying continuum Schwarzschild 
space-time and its geodesics therefore become more sim¬ 
ilar to those in the continuum space-time. 

Thus combined with our analytic arguments provided 
at the end of Section IVTl we conclude that (1551) gives the 
correct velocity as a function of r for particles following 
radial time-like geodesics. 

APPENDIX C: RADIAL LENGTHS OF 
SCHWARZSCHILD REGGE BLOCKS 

In our simulations of the flat and open universes, 
our Regge blocks’ radial lengths were not constant but 
were increased as the blocks were further away from 
the Schwarschild-cell centre. This decision was moti¬ 
vated by the fact that the underlying space-time being 
approximated becomes increasingly flat. We therefore 
attempted to maintain higher resolution while the cur¬ 
vature was high but then decrease the resolution with 
minimal impact on accuracy as the curvature decreased. 
This technique allowed us to perform larger-scale sim¬ 
ulations within more attainable computation times. We 
now describe our method for specifying the block’s radial 
length. 

In Schwarzschild space-time, the radial distance be¬ 
tween radius r'j > 2m and some arbitrary r > r% is given 
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by 



L J x = Vi 

(Cl) 

When r = 1 , this is identical to (l5lT) . Distances along 

‘radial edges’ of Regge blocks approximate radial dis¬ 
tances in Schwarzschild space-time by linearly interpo¬ 
lating d(r) between ri and r,; + i, as shown in Fig. [35] 
Since d(r) is convex, the interpolation will always under¬ 
estimate the true distance, and the error e of the approx¬ 
imation is given by e = max [d(r) — lin(r)], where lin(r) 
denotes the interpolating function. 



FIG. 35. The Regge block interpolation intersects the 
Schwarzschild graph at r = n and r = n+i, where ri < n+ i 
and where ri corresponds to the origin of the graph; it 
interpolates the radial distance for only those values of r 
lying between these two points of intersection. Since the 
Schwarzschild distance function is convex, the interpolation 
will always underestimate the distance in this regime. 

Our goal is as follows. Given some error-tolerance £toi, 
we want our blocks to be as long as possible in the radial 
direction while still satisfying e < E to i- If we enter a 
new block from its left or right face, then we would need 
to re-calculate the block’s radial length h = rj+i — r, 
in accordance with our goal. It will always be the case 
that we know one of r,; or r"j +\ depending on which face 
we entered by, and we would need to deduce the other 
quantity. However our computation must give consistent 
results for h regardless of which quantity we used. In 
our particular implementation, we have also chosen to 


impose the constraint that h must be some integer units 
of a minimal interval Aro, a parameter which we can 
freely specify. 

The most natural approach would be to solve d'(r) = 
lin^r), which maximises d{r) — lin(r), to obtain r = r max , 
the point at which the error is greatest. Also through its 
dependence on lin(r), r m ax would be a function of h. We 
would then set e = e to i and solve d(r max ) - lin (r ma x) = 
£ to i for h. However, this second equation can only be 
solved numerically, and from a programming point of 
view, it was easier to implement a different approach in¬ 
stead. 

We note that for r > r^, the functions d{r) and lin(r) 
are bounded from above by 

d + (r) = [( r - A) + d(ri), (C2) 

since this is just the equation of the tangent to d(r) at 
ri. And also, for r > r*, the two functions are bounded 
from below by 

d~(r) = (r - n) + d(ri); (C3) 

this follows because the gradient of d(r) goes asymptot¬ 
ically to unity from above as r —> oo, which means d{r) 
will always increase more quickly than d~ (r). Therefore, 
after the two functions have intercepted at r — n, d{r) 
will always be strictly greater than d~(r). The func¬ 
tion lin(r) also intercepts both d(r ) and d~(r ) at r = r^. 
But it intercepts d(r) again at r.; + i > ri, a point where 
d(r) > d~ (r). Because lin(r) is a linear function, this in¬ 
terception at r ^|-i implies that lin(r) must have a greater 
gradient than d~(r), and therefore lin(r) must also be 
strictly greater than d~(r) in the region r > r;. 

Thus we can bound our error by 

e < max [d + (r) - d~(r)] = d + (r i+1 ) - d~(r i+1 ). 

We shall therefore require that 

£toi > d + {r i+ 1 ) - d~(r i+ i), 

and solve for r,_|_i from this equation; this yields the so¬ 
lution 

7 _ £t0l 

\/Til (ri - 1 ) - 1 

To get ft as a function of r,; + ] instead, we set r,; to be 
n = ri +1 — h and solve for h in the preceding equation; 
this yields 


h = 


&tol 


2 ( 2 e tol - 1 ) 


2?y+i + Etoi — 2 + 


4r i+ i(r i+ i 


1 — Etol) + £toi(4 + Etol) 
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As mentioned above, we want h to satisfy h = n Aro 
for some integer n. Therefore our actual h is obtained by 
taking n to be 


£tol 


Ar 0 (\Z r i/{n - 1) - l) 


(C4) 


or 


£tol 


_2Ar 0 {2e to i - 1) 


2r^+i + £toi — 2 

4r i+ i(r i+ i - 1 - Etoi) (C5) 
11/2 > 

+ £toi( 4 + £tol ) 


as appropriate, where [x\ gives the greatest integer less 
than or equal to x, and then substituting this n into 
h = n Aro. 

Finally, we want to ensure that our algorithm would 
give the same h regardless of whether it used r* or rj+i. 
To satisfy this requirement, whenever the program cal¬ 
culated Vi+\ from ri, it would check to see if it could 
recover the same h using the new value for n+ 1 - If not, 
then it would decrement h by one unit of Aro and check 
again, and it would continue decrementing until obtain¬ 
ing agreement. The program can only decrement if the 
error is to remain within tolerance. We had also required 
that h be at least Aro; thus the decrementing would stop 
if h reached this length. 


APPENDIX D: BOUNDARY CONDITIONS 
FOR SCHWARZSCHILD-CELLS OF 
UNEQUAL MASSES 

An observer sitting at the interface of two cell bound¬ 
aries must observe the same physics regardless of which 
co-ordinate system the observer uses. In particular, the 
following two conditions must be satisfied: 

1 ) the observer must measure the same proper time 
regardless of which cell’s metric is used; 

2 ) and the observer must measure the same spatial 
distances locally along the boundary regardless of 
which cell’s metric is used. 

However if the cells in the lattice are no longer identical, 
we would no longer have the same lattice symmetries as 
before. In particular, the exact location of the boundary 
between two unequal cells becomes less transparent, as 
it is no longer equidistant to the two cell centres. Re¬ 
call that in the perfectly symmetric lattice, a test parti¬ 
cle sitting on the boundary between two cells would by 
symmetry always remain at the boundary and yet fall 
simultaneously towards both cell centres. Although the 
symmetry is no longer present when the two cells are no 


longer identical, we shall still assume that test particles 
at the boundary fall simultaneously to both centres. This 
condition defines the location of the cell boundary for us: 
it is effectively defined to be where the two cells’ gravita¬ 
tional influences are equal. On one side of the boundary, 
the gravitational influence of one mass dominates, and 
we approximate the space-time by a Schwarzschild space- 
time centred on that mass. On the other side, the other 
mass dominates, and we approximate the space-time with 
a Schwarzschild space-time centred on that mass. Once 
again, we can understand this simultaneous free-fall of 
the boundary towards the two centres as actually the 
motion of the two masses towards each other under their 
mutual attraction. This mutual attraction gives rise to 
the expansion and contraction of the lattice itself, which 
manifests as the expansion and contraction of the cell 
boundary. Given this definition and our two boundary 
conditions, we can now derive a set of constraints that 
E b and the cell radius r b must satisfy. 

Local spatial distances along the boundary are given 
by 


r b dtt. (Dl) 

Since by assumption, particles co-moving with the 
boundary follow radial geodesics, then this distance 
evolves as 


A dsl . 

dr 


By condition (2) we require that 

r bl dfli = r b2 dVL 2 , 


(D2) 


(D3) 


where the numerical subscripts refer to cells 1 and 2; 
combined with condition <H3 we additionally require that 


dr dr 


(D4) 


since dr\ = dr 2 = dr. Combining these two conditions, 
we obtain 


1 dr bl 1 dr b2 


r bl dr r b2 dr 

and using (Hhd, we can express this equivalently as 


IE bl - 1 + 


2 mi 


I E b2 ~ 1 + 


2m 2 

r b2 


(D5) 


(D6) 


By condition (1) we require that Ati = At 2 . Recall 
that the proper time of a freely falling particle is given 
by equations EU1) to C3l) . However, from the form of 
equations to (1221) . it is clear that we cannot satisfy 
this condition unless both cells are of the same type, that 
is, both cells are open, flat, or closed. This constrains 
the cases we need to consider to just three, and we shall 
consider each in turn. 
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If E b = 1 for both cells, then (l2lT) and (l23l) imply that 


Ari = A r 2 


2 1 ^ 3/2 

3 


2 1 
3 y/2m 2 



- T 0 , 


where To is a constant of integration. From this, we ob¬ 
tain the relation 


n 2 = 



(D7) 


It can be checked that this relation also satisfies (lD6l) and 
hence condition (2) as well. 

If E b < 1 for both cells, then (120[) and (l23l) imply that 


2toi 


(l-£ bl ) 3/2 


1 ~E bl 


2toi 
2 m 2 


r bl \ 1 - 


1 ~ E bl 

2rai 


r bl + cos 


-1 . / I -E bl 
2mi 


r bl 


- t 0 


(1 ~E b2 f 2 


1 ~ E b2 
2m 2 


n 2 \ ! - 


1 — E b 2 
2m 2 


r b2 + cos 


1 ~ E b 2 
2m 2 


-n 2 


to- 


I- 

The two sides of this equation can be made equal if we The second equation leads to the constraint 
simultaneously equated 


1 ~ E bl 

2 toi 


1 — Ei 


-Tb! = 


b2 


2 m 2 


-n 2 


and 


2mi 


2 m 2 


(1 ^E bl f 2 (l~E b2 ) 3/2 ' 


' \ 2/3 

(E b2 (E bl - 1), 


(D8) 


and if we substitute this into the first equation, we re¬ 
cover (lD7l) . Again, it can be checked that (1D71) and (ID8I) 
combined satisfy m and hence condition (12)] as well. 
Finally if E b > 1, then (l22ll and (l23l) imply that 


2mi 


(E bl - 1) 3/2 


Et!-! 


2rrii 
2TO2 


r bi \ 1 + 


E bl - 1 

2toi 


(- E b2 - 1) 3/2 


E b2 - 1 


2 m 2 


r b 2 \ ! + 


r b , — sinh 


E b2 - 1 
2 m 2 


:„U-1 Eb!-!. 


2mi 


-r bl 


- T 0 


r b2 — sinh 


-l 


E b2 - 1 
2 m 2 


-n 2 


- To ■ 


As with the E b < 1 case, the two sides of this equation 
can be made equal if we simultaneously equated 

E bl — 1 E b2 — 1 

—T- r bl = — - r b2 

2toi 2 m 2 


and 


2mi 2 m 2 

(E bl -1) 3/2 “ (E b2 — 1) 3/2 ' 

This clearly leads to the same constraints as in the E b < 1 
case. 

We therefore find that for all cases, two neighbour¬ 
ing cells must satisfy constraints (ID7D and (ID8I) at the 
boundary^ And when m i = m 2 , we recover E bl = E b2 
and r bl = r &2 . 


18 The author wishes to acknowledge that constraints llDTll and 
(1D8I) were actually first derived by Ruth Williams in an unpub¬ 
lished calculation, where she derived the boundary conditions us¬ 


ing the Israeli junction conditions instead. Equating the two met¬ 
rics on the boundary gave the condition dr^/dri — drb 2 /dr 2 , 
and equating Tr where k is the extrinsic curvature, gave con¬ 
dition lD5l l. 
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